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ABSTRACT 

We have derived isotopic fractions of europium, samarium, and neodymium in two metal-poor giants 
with differing neutron-capture nucleosynthetic histories. These isotopic fractions were measured from 
new high resolution (R ~ 120,000), high signal-to-noise (S/N~ 160-1000) spectra obtained with the 
2dCoude spectrograph of McDonald Observatory's 2.7 m Smith telescope. Synthetic spectra were gen- 
erated using recent high-precision laboratory measurements of hyperfine and isotopic subcomponents 
of several transitions of these elements and matched quantitatively to the observed spectra. We inter- 
pret our isotopic fractions by the nucleosynthesis predictions of the stellar model, which reproduces 
s-process nucleosynthesis from the physical conditions expected in low-mass, thermally-pulsing stars 
on the AGB, and the classical method, which approximates s-process nucleosynthesis by a steady 
neutron flux impinging upon Fe-peak seed nuclei. Our Eu isotopic fraction in HD 175305 is consistent 
with an r-process origin by the classical method and is consistent with either an r- or an s-process ori- 
gin by the stellar model. Our Sm isotopic fraction in HD 175305 suggests a predominantly r-process 
origin, and our Sm isotopic fraction in HD 196944 is consistent with an s-process origin. The Nd 
isotopic fractions, while consistent with either r-process or s-process origins, have very little ability 
to distinguish between any physical values for the isotopic fraction in either star. This study for 
the first time extends the n-capture origin of multiple rare earths in metal-poor stars from elemental 
abundances to the isotopic level, strengthening the r-process interpretation for HD 175305 and the 
s-process interpretation for HD 196944. 

Subject headings: atomic data — nuclear reactions, nucleosynthesis, abundances — stars: abundances 
— stars: Population II 



1. INTRODUCTION 

A wealth of observational studies over the last decade 
have focused on measuring precise abundances of neu- 
tron (n)-capture elements in metal-poor stars. The heart 
of this enterprise is identifying the origin of the nuclides 
with Z > 30 in our Galaxy. These nuclides are produced 
in stars through either the rapid (r)-process or the slow 
(s)-process. Rapid-process nucleosynthesis occurs in en- 
vironments with high neutron fluxes and densities, al- 
lowing many neutrons to be captured by existing nuclei 
much more rapidly than the timescales for /3-decay. Slow- 
process nucleosynthesis occurs where the neutron densi- 
ties are lower, with /3-decay timescales generally shorter 
than the times between n-capture events. To build heavy 
nuclei through the s-process, a chain of stable or long- 
lived nuclei must exist from the existing seed nucleus to 
the final product. The r-process and s-process lead to the 
creation of different sets of heavy nuclei, some of which 
can only be created in one process or the other ("pure- 
s-" or "pure-r-nuclei" ) and some of which are created by 
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both processes (see, e.g.. lCowan fc Snededl2006T) . While 
the precise astrophysical site(s) of the r-process have yet 
to be identified, s-process nucleosynthesis takes place in 
low- and intermediat e-mass stars on the asymptotic giant 
branch (AGB) (e.g., [Eisoet^Q]!!). 

To illustrate the effects of both n-capture processes, 
in Figure [T] we show a table of nuclides covering por- 
tions of the atomic number domain 60 < Z < 63. The 
s-process can add only one neutron before a /3-decay oc- 
curs for an unstable nucleus. On the other hand, the 
r-process quickly floods a target nucleus with many neu- 
trons, overwhelming the /3-decay rates. When neutrons 
are no longer being added, the nucleus will /3-decay re- 
peatedly until it reaches stability. For isotopes produced 
primarily in the s-process (e.g., 146 Nd, 148,150 Sm) a sta- 
ble isotope with one less nuclide always exists. For iso- 
topes produced primarily or wholly by the r-process (e.g., 
150 Nd, 147,149,152, i54 Sm ) ; there ig & clear path for ^ecay 

from unstable nuclei to the lower right (in the table of 
nuclides) of the stable r-process isotope, sometimes un- 
reachable by the s-process. 

The nucleosynthetic signatures of the s- and r- 
processes stand in sharp contrast to one-another. The 
Solar System (S. S.) isotopic abundances are well- 
established fro m studies of CI m e teorite s, as summa- 
rized in, e.g., lAnders fc Grevessel |l989) and lLoddersI 
(2003]). Two' primary methods are used to decompose 
the S. S. isotopic abundances into their constituent s- 
and r-process origins. In the "standard" o r "clas sical" 
method, firs t artic ulated by IClavton et al.l (|1961h and 
ISeeger et all (119651) and w ith abundances later u pdate d 
by. e.g.. iKappeler et all (|1989l ): iBurris et all f|2000f >: 
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ISimmerer etall (|2004l ). and ICowan et all (|2006l) . the 
s-process is modeled as a slowly- varying function of neu- 
tron exposure. This method attempts to smoothly fit 
the u aN s " curve (the product of the n-capture cross- 
section and the S. S. s-process abundances) using lab 
measure ments of a to determ ine N s . In the "stellar" 
model of lArlandini et al.1 (|1999h . isotopic abundances for 
s-process nucleosynthesis in 1.5— 3.0 M Q stars are pre- 
dicted based upon nuclear physics cross sections and stel- 
lar model calculations for thermally-pulsing (TP) stars 
on the AGB. In both of these cases, the r-process abun- 
dances are then derived as residuals between the total 
S. S. abundances and the s-process products. Both ap- 
proaches have strengths and weaknesses. The classical 
method is model-independent, but it assumes an em- 
pirical, smoothly-varying relationship for the <jN s abun- 
dance curve (which relies on the measured S. S. abun- 
dances) and ignores details of the nuclear physics. The 
stellar model does not depend on knowledge of the S. S. 
abundance distribution, but is heavily dependent on nu- 
clear physics laboratory measurements and complex stel- 
lar AGB model calculations. 

These models can also be used to predict the rela- 
tive amounts of s- or r-process nucleosynthesis prod- 
ucts in metal-poor stars. Several key rare earth ele- 
ments are commonly used to assess r-/s-process abun- 
dance dominance in stars. 6 Eu is commonly used 
as a surrogate for r-process nucleosynthesis because 
it is so predominantly composed of r-process material 
and has several easily-observed spectral lines. Simi- 
larly, Ba or La are commonly used as surrogates of 
s-process nucleosynthesis. Thus, many studies have em- 
ployed [Ba/Eu] 7 ratios in metal-poor stars to estimate 
the re lative contribu t ions f ro m the s- and r-pr o cesses 
(e.g.. ISpite fc Spitd [1975 iGratton fc Snedenl [l99l 
Burris et al. 2000; Mashonkina et al. 2003; Honda et al. 
20041* iBarklem et all [2001 Francois 1TaT]l2007[ L which 
in turn const rain models of chem i cal evolution in the 
early Galaxy (iTrayaglio et al.lll999l l2001h ICescutti et "all 
120061 iKratz et al.ll2007f l. Recently. [La/Eu] ratios have 
been employed as an alternative to [Ba/Eu] to avoid 
difficulties in Ba abundance determinations (see, e.g., 
Simmerer et~aT1l2004t Uonsell et all2006l : IWinckler et all 
III). 

Elemental abundances, the sums of abundances of their 
naturally-occurring constituent isotopes, are more easily 
measured in stars than the isotopic abundances; however, 
isotopic abundances should be more fundamental indica- 
tors of n-capture nucleosynthesis because they can di- 
rectly confront r-process and s-process predictions with- 
out the smearing effect of multiple isotopes. Much recent 
evidence (e.g., the growing number of r-process-enhanced 
stars that conform to the scaled-S. S. r-process elemen- 
tal abundance distribution) supports the hypothesis of a 
universal r-process mechanism for elements with Z > 56. 

6 The s-process percentages for some of the commonly-studied 
rare earth elements, predicted by the stellar model and classical 
method, respectively, are: Ba (five naturally-occurring isotopes), 
81.1% or 85.3%; La (one isotope), 62.1% or 75.4%; Nd (seven iso- 
topes), 58.8% or 57.9%; Sm (seven isotopes), 30.7% or 33.1%; Eu 
(two isotopes), 5.8% or 2.3%. 

7 We adopt the usual spectroscopic notations that [A/B] = 
logi (N A /N B )* - logio(N A /N B )Q, and that log e(A) = 
logi (N A /N H ) + 12.00, for elements A and B. 



It is important to demonstrate that the r-process elemen- 
tal abundance pattern extends to the isotopic level. Such 
agreement would greatly strengthen the argument for a 
universal r-process mechanism for the heavy n-capture 
elements. 

Most isotope fractions, unlike elemental abundances, 
are very insensitive to the model atmosphere parameters. 
The wavelength of a spectral line is split by two effects, 
hyperfine structure (HFS) and isotope shifts, and the 
isotopic fractions can be measured by detailed compar- 
isons of an observed absorption line profile to synthetic 
spectra of these line substructures. If this splitting is 
comparable to or greater than the combined effects of 
stellar thermal, macroturbulent, and rotational broad- 
ening and the broadening caused by the spectrograph 
slit, then it may be possible to measure the isotopic frac- 
tions. In principle, unlike with elemental abundances, 
any n-capture element with multiple naturally-occurring 
isotopes that are produced in different amounts by the 
s- and r-processes can be used to assess the relative s- 
and r-process contributions to the stellar composition. In 
practice, the availability of quality atomic data has lim- 
ited studies of isotopic fractions of rare earth elements 
in stars, but the relatively recent increase in laboratory 
spectral line substructure studies has led to the identi- 
fication of a few lines of a few elements that might be 
a nalyzed at the is o topic level. 

iCowlev fc Frevl (fl989h and iMaeain fc Zhaol (|1993ah 
first suggested that the width of the Ba II 4554A line 
in stellar spectra would be sen si tive to t he is otopic 
fraction. iMagain fc Zhaol (ll993aT) . IMaeain! (|1995h . and 
lLambert fc Allende Prietol (|2002f ) all examined the Ba 
isotopic fraction in the metal-poor subgiant HD 140283. 
Their measured isotopic fractions had uncertainties too 
large to unambiguously discern the Ba nucleosynthetic 
origin in this star. A rec ent series of papers by 
Mash onkina et all (|1999t 120031 ) and IMashonkina fe Zhaol 
(2006) have measured Ba isotopic fractions in thin and 
thick disk stars. Their approach is to set the Ba el- 
emental abundance from lines not affected by HFS or 
isotope shifts and then to derive the isotopic fraction 
from the 4554A line, whose shape is sensitive to both 
its elemental abunda n ce an d isotopic composition of Ba. 
iMashonkina fc Zhaol ([2006) clearly show that a higher 
fraction of r-process Ba isotopes is present in thick 
disk stars than in S. S. material and thin disk stars, 
and they also show that a smooth trend of increasing 
s-process contributions occurs with increasing metallic- 
ity Since only the 4554A line is useful for examining 
the Ba isotopic mix, and the profile of this line is known 
to suffer from non-LTE e ffects in metal-poor stars (e.g., 
IMashonkina et all 119991 IShort fc HauschildtJ |2006h . we 
do not attempt to assess the Ba isotopic fractions in this 
study. The Mashonkina et al. studies of the Ba isotopic 
fraction do, however, employ non-LTE calculations. 

Eu isot opic fractions in met al-poor stars were first re- 
ported bv ISneden et al.l (|2002h . who found that the iso- 
topic ratio in three r-process-rich giants was consistent 
with the S. S. iso topic ratio. Subsequent studies by 
lAoki et ail (|2003al lbl) found that Eu isotopic ratios could 
marginally distinguish s- and r-process contributions in 
stars that had prev iously been shown to exhibit such 
chemical signatures. iLundavist et alj (|2007f ) have inves- 
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tigated the Sm isotopic fractions in a very metal-poor 
r-process-enriched star, qualitatively showing that the 
isotopic mix is in agreement with an r-process origin. 
No measurements of the Nd isotopic fraction have been 
made outside of the S. S., and no study has yet attempted 
to measure isotopic fractions of multiple n-capture ele- 
ments in the same star. 

In our study, we analyze the isotopic fractions of Eu, 
Sm, and Nd in two metal-poor giant stars. Our obser- 
vations and methods of analysis are described in § [5] and 
IH respectively. We summarize the available atomic data 
for each species in § |3l Our measurements of the ele- 
mental abundances and isotopic fractions are presented 
in § [5] and §0 In § we discuss the implications of 
our measurements in the context of using isotopic frac- 
tions of multiple species to constrain the s- and r-process 
nucleosynthetic histories. In § [5] we use Eu and Pb 
abundances collected from the literature to assess the 
s-process model predictions for the 151 Eu isotope, since 
pure-r- and pure- s-process nucleosynthesis reactions pro- 
duce such contrasting abundances of these elements. We 
also make suggestions for how future studies might uti- 
lize measurements of n-capture isotopic fractions to make 
similar comparisons for other isotopes in § [HI We sum- 
marize our results in § 1101 In the appendices, we in- 
clude comments describing the nuances associated with 
analyzing each line of each species (Appendix |A|). and 
we include our computations of the hyperfine structure 
subcomponents for both Sm (Appendix IB]) and Nd (Ap- 
pendix [C|). 

2. OBSERVATIONS 
2.1. Selection of Target Stars 

We choose two bright, metal-poor giants for our anal- 
ysis, HD 175305 and HD 196944. The technical require- 
ments of achieving high signal-to- noise (S/N) ratios at 
very high resolution on a telescope with a moderate aper- 
ture mandate that we select bright (V < 8.5) targets. We 
select metal-poor giant stars for this analysis because 
(1) metal-poor stars present significantly less crowding 
of strong absorption features in the blue spectral region 
where the lines of interest are concentrated, (2) absorp- 
tion lines of all species will exhibit greater line depth in 
giants than in dwarfs because of the decrease in the con- 
tinuous opacity with decreasing electron pressure in gi- 
ants, and (3) at a given metallicity, molecular formation 
will be less in giants than in dwarfs. To maximize the 
expected contrast in the isotopic fractions, we have delib- 
erately selected stars with different chemical enrichment 
histories. Although these stars have been well-studied, 
no previous attempts have been made to examine their 
isotopic fractions. 

HD 175305 is a metal-poor ([Fe/H] = -1.6), bright 
(V ~ 7.2) giant that was selected because it shows 
evidence of r-process enrichment ([Eu/Fel ~ +0.5, 
[Ba/Eul < 0.0. [La/Eu] ~ -0.4 ; iBurris et all 120001 : 
lFulbrighdl200(l iCowan et all 12001 . HP 175305 is not 
enric hed in carbon ([C/Fe] = +0.2; IWallerstein et alJ 
1979). Although only a handful of heavy n-capture 
abundances have been measured previously, the high 
[Eu/Fe] and low [Ba/Eu] ratios are suggestiv e of an 
r-process origin. Combining 1 Cowan et al.l (|2005f )'s abun- 
dance measurements of osmium (Os, Z = 76), iridium 
(Ir, Z = 77), and platinum (Pt, Z = 78) from UV spec- 



tra and lLawler et all ((2007) "s measurement of hafnium 
(Hf, Z = 72) with additional heavy n-capture elemental 
abundances measured in our spectra will produce a more 
thorough understanding of this star's enrichment history 
to complement our isotopic fraction measurements. 

HD 196944 is a very metal-poor ([Fe/H] = -2.4), 
bright (V ~ 8.4) giant with clear s-process signatures. 
This star has high C, N, and O abun dances ([C/Fe] ~ 
+1.3, [N/Fe] ~ +1.3, [O/Fe] ~ +1 4; iZacs et all 1199^ 
lAoki et all 120021 ; Uonsell et ai1l2005t ). high Ba and La 
abundances ([Ba/Fe] ~ +1.1, [La/Fe] ~ +0.9), a low 
Eu abundance ([Eu/Fe] ~ +0.2) and a large over- 
abundance of lea d f[Pb/Fe] ~ +2.0:lVan Eck et aLl l200lL 
lAoki et al.| [200l IVan Eck et alj 12003ft. HD 19694 4 also 



exhibits radial velocity variations ( Lucatello et alJ r2005) , 
indicating that it is a member of a binary or multiple star 
system. These characteristics have led to the classifica- 
tion of HD 196944 as a carbon-enha nced metal-poor s- 
enrich e d (CEMP-s) star, as d efined in lBeers fc Christliebl 
(|2005ft . iGallino et all (120051) have shown that the abun- 
dance pattern of HD 196944 can be reproduced by as- 
suming mass transfer of material enriched in s-process 
elements from a companion star with an initial mass of 
1.5 M© that passed through the TP-AGB phase. 

We also observed HD 122563 in our study; at V ~ 6.2, 
this star is the brightest very metal-poor ([Fe/H] = —2.8) 
star in the sky. Unfortunately the lines of Eu, Sm, and 
Nd proved to be too weak to provide any assessment 
of the isotopic fractions. This is unfortunate because 



Honda et all (120071), following studies of this star by 
Cowan et all (|2005l) . lAoki et all (|2005h . and lHonda et all 
fl2006j), have suggested that HD 122563 and HD 88609 



may be the first two known representatives of a class of 
stars whose n-capture enrichment has been dominated 
by the weak r-process, which produces a distinctly dif- 
ferent abundance pattern than the main r-process. 8 Our 
analysis of HD 122563 did allow us to assess the location 
and strength of blending features for some of our lines of 
interest. It is possible that the Eu n lines at 3819, 3907, 
and 4129A, which were not covered in our spectrum, may 
give some hint to the Eu isotopic fraction. 

2.2. Description of Observations 

For a set of lines of a given species with the same ion- 
ization state in a stellar atmosphere, the relative strength 
of each line depends on its transition probability, excita- 
tion potential, wavelength, and the continuous opacity of 
the atmosphere at that wavelength (see, e.g.. lGravll2005l . 
p. 389). In metal-poor stars the continuous opacity (from 
H — ) increases slowly over the visible wavelength region as 
the wavelength increases; these two terms take opposite 
signs in calculations of the line strength and roughly off- 
set each other. To first approximation then, the relative 
strengths of individual absorption lines are determined 
by the t ransition properties of the lines themselves. Fig- 
ure 10 of lLawler et al.l (|2Q06h illustrates this point for the 
case of Sm 11 lines. For the n-capture species considered 
in our study, the first-ionized state is dominant in our 
stars, and most of the transitions have similar excita- 
tion potentials (~ 0.0-0.7 eV), so the oscillator strengths 



8 P lea se see, e.g., [W asscrbur g et all 11999. IQ ian fc Wasserbur 
200(1 IWasserburg fc Qianl 



2000, 



Johnson & Boltc 200 



Wanaio fc Ishimarul 120061 . and IKratz et al.1 [2007 for further 



discussion of these two processes. 



4 



Roederer et al. 



most significantly affect the relative line strength. The 
log(g/) values generally decrease when moving from the 
near-UV to the near-IR spectral regions, so the strongest 
lines in stellar spectra tend to lie in the near-UV and 
blue end of the visible spectrum. This is unfortunately 
where the most severe blending with other spectral fea- 
tures is likely to occur. The isotope shifts increase with 
increasing wavelength; this broadens the stellar line pro- 
files, making it easier to measure the isotopic fraction. 
The decrease in line strength, decrease in blending fea- 
tures, and the increase in isotopic shifts at longer wave- 
lengths compete against each other, producing a range in 
the visible spectrum from ~ 4400— 5100A where isotopic 
fractions are most easily measured. 

We acquired new observations of our target stars us- 
ing the 2dCoud e cross-dispersed echelle spectrograph 
(|Tull et al.lll995D on the 2.7 m Harlan J. Smith Telescope 
at the W. J. McDonald Observatory. Observations of 
HD 122563 and HD 196944 were made on 2006 June 12- 
18, and observations of HD 175305 were made on 2006 
August 08-12. Sets of 30- minute exposures were taken 
with a Tektronix 2048 x 2048 CCD ("TK3"). Around 
5000A, to sample the full width of a line profile (~ 0.2A) 
with at least 10 points requires AA ~ 0.02 per pixel and 
a two-pixel resolution of R ~ 125,000. This requirement 
dictates the spectral resolution necessary to complete 
our analysis. Our data have FWHM resolving powers 
R ~ 120,000. Lines in our Th-Ar comparison spectra 
can be fit with a pure Gaussian function, and no variation 
in the instrumental profile was found from night-to-night 
over the course of either observing run. Due to the lim- 
ited wavelength coverage at this setup, we have observed 
each star with three grating settings to increase the wave- 
length coverage. Our spectra cover 4200 < A < 6640A 
in HD 122563 and HD 196944 and 4120 < A < 5890A 
in HD 175305; however, even when using multiple grat- 
ing settings our wavelength coverage is still incomplete 
within these ranges. 

Image processing, order extraction, and wavelength 
calibration were performed using standard IRAF 9 tasks 
in the echelle and onedspec packages. Echelle orders were 
cross-correlated and co-added, and the continuum nor- 
malizatio n was performed with the SPECTRE software 
package (|Fitzpatrick fc Snedenlll987f ). Our final S/N val- 
ues (per pixel) range from ~ 160 at 4130A to - 330 at 
5100A in HD 175305, - 160 at 4200A to - 400 at 5000A 
in HD 196944, and ~ 430 at 4200A to ~ 1000 at 5000A 
in HD 122563. 

3. ATOMIC DATA 

3.1. HFS Structure and Isotope Shifts 

HFS is primarily caused by the magnetic interaction 
between an electron and the nucleus of an atom. This in- 
teraction occurs when the nucleus has a magnetic dipolc 
moment that results from a non-zero angular momentum, 
which occurs when the nucleus has an odd number of pro- 
tons and/or an odd number of neutrons. The strength of 
this interaction is characterized by the magnetic dipole 
"A" constant, which can be measured in laboratory stud- 

9 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



ies. We also include the measured electric quadrupole 
interaction "B" constant in our HFS computations, but 
magnetic interaction is usually dominant. In stellar spec- 
tra, HFS is manifest as a broadening of a single line of a 
particular atom due to the energy shifts. Isotope shifts, 
on the other hand, are slight rearrangements of the en- 
ergy levels of one isotope of a particular species relative 
to another isotope of the same species. The most use- 
ful lines of heavy atoms — such as those considered in this 
study — usually have isotope shifts dominated by the field 
shift. This shift is from the finite volume occupied by the 
nucleus of an atom. Proton charge is distributed within 
the nucleus, and not concentrated at the origin, hence 
the electric field inside the nucleus does not have a 1/r 
dependence. For electrons that have a significant prob- 
ability of existing at small r, such as the 6s electron in 
the ground state configuration of Eu II, this can have a 
sizable effect. Both the reduced mass shifts and specific 
mass shifts of the upper and lower levels also contribute 
to the total isotope shifts of a line. Isotope shifts for our 
elements have been measured in laboratory experiments. 
In stellar spectra, the effect of isotope shifts is to broaden 
the absorption line profile when more than one isotope 
of the same atom are present in the stellar atmosphere. 

3.2. Europium Atomic Data 

There are two naturally-occurring isotopes of Eu: 
151 Eu (47.8% S. S.) and 153 Eu (52.2%). Both Eu iso- 
topes exhibit HFS. Both isotopes of Eu are produced in 
roughly equal amounts in the r-process, which dominates 
over the s-process as the origin of the S. S. Eu. Because 
there are only two Eu isotopes, it is straightforward to 
measure the isotopic fraction by measuring the fraction 
of one of these isotopes relative to the total amount of 
Eu. The fraction /151 is defined as 



/151 — 



iV( 151 Eu) 



(1) 



AT(Eu) 

In Table [2j we list the fraction of 151 Eu that would be 
expected if all the Eu present in a given star were to 
have originated in only the s-process or the r-process for 
both the stellar model and classical method. The stellar 
model predicts a /* 51 abu ndance that is consistent with 
two (r + s)-enriched stars (|Aoki et al.ll2003bl 120061) . but 
the classical method predicts a /f 51 abundance that is 
nearly 30 times smaller. We also note that even the pre- 
dictions for the pure-s-process abundance of 151 Eu vary 
greatl y within different sets of classical method calcula- 
tions (jArlandini et al.l 119991 : ICowan et ani2006h . so this 
discrepancy is not limited to only the classical method 
versus stellar model predictions. Lengthy discussion of 

this discrepancy will be taken up further i n § [8] 

We a dopt the Eu 11 log(gf) values of lLawler et alJ 
(2001b]). The hyperfin e and isotopic components are 
taken from llvans et all (|2006l) . The HFS of the three 
Eu 11 lines that we consider is very wide, ~ 0.2— 0.3 A, 
which facilitates our ability to measure the isotopic frac- 
tion. 

3.3. Samarium Atomic Data 

There are seven naturally-occurring isotopes of Sm: 
144 Sm (3.1% S. S.), 147 Sm (15.0%), 148 Sm (11.2%), 
i49q™ mq so/a i50 Sm ( 7 .4%) ; i52 Sm ( 2 6.7%), and 154 Sm 



Sm (13. 
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(22.8%). The 144 Sm isotope is exclusively produced by 
the p-process; fortunately for this analysis, this isotope 
only comprises 3.1% of the Sm in the S. S., so we can 
neglect its contribution to the isotopic mix. According 
to the stellar model and classical method, 148 Sm and 
Sm are produced exclusively by the s-process, while 
the r-process dominates the production of the other four 
isotopes. 

The two odd- A isotopes exhibit HFS structure in addi- 
tion to the energy shifts between all seven isotopes. The 
Sm isotope shifts are ~ 0.05 — 0.10 A, which is slightly 
larger than the Sm HFS splittings, ~ 0.05 A, yet both 
the Sm isotope shifts and HFS splittings are significantly 
smaller than the Eu HFS splittings. The two heaviest 
isotopes, 152 Sm, and 15 Sm, have large isotope shifts rel- 
ative to the isotope shifts of the lighter five isotopes in 
ma ny of the transitions we analyze here (cf. Figures 1 and 
2 of lLundqvist et al1l2007t ). Therefore, it makes sense to 
define the isotopic fraction in terms of these two heaviest 
isotopes. The quantity /152+154 is denned by 



152+154 



N( 152 Sm) + 7V( 154 Sm) 
AT(Sm) 



(2) 



In Table [2 we also list the isotopic fraction of Sm that 
would be expected if all the Sm present in a particular 
star were to have originated in only the s-process or the 
r-process for both the stellar model and classical method. 
Note that this definition more easily distinguishes pure-s 
and pure-r isoto pes of Sm than the /odd definition sug- 
gested in § 7 of lLawler et al.l ((2006). In addition, the 
pure-s and pure-r isotope fractions predicted by the stel- 
lar model and classical method do not differ significantly 
when expressed in terms of /i52+i54- 

We have computed the HFS structure patterns for 13 
Sm lines in this wavelength range usin g the hyperfine A 
and B constants and isotope shifts of iMasterman et all 
(2003), whe n available, supplemen ted with additional 
values from iLundqvist et al.l (|2007f) and others as de- 
scribed in Appendix IB] Three of these lines (4719, 5052, 
and 5103) do not have complete sets of hyperfine A and B 
constants reported in the literature. To gauge the impact 
of the HFS splitting on our stellar syntheses, we have 
compared two syntheses of lines with full sets of hyperfine 
A and B constants, one including the HFS splitting and 
isotope shifts and one including only the isotope shifts. 
Syntheses with a lower fraction of /152+154 are slightly 
more broadened in the case of full HFS treatment than 
when only the isotope shifts are considered. The differ- 
ences are small; the synthesis including HFS is f» 5—6% 
wider at FWHM than the synthesis with no HFS (when 
the percentage of isotopes exhibiting HFS is maximized 
in the synthesis). Nevertheless we recommend that lines 
lacking complete sets of hyperfine A and B constants 
should not be used to measure the isotopic fraction by 
fitting the shape of the line profile; however, they should 
still be useful for other methods of measuring the iso- 
topic fraction and for elemental abundance analyses and 
are included in the discussion he re for this reason. We 
adopt the Sm 11 log(gf) values of lLawler et alj (|2006). 

3.4. Neodymium Atomic Data 

Nd, like Sm, has seven naturally-occurring isotopes in 
S. S. material: 142 Nd (27.2% S. S.), 143 Nd (12.2%), 144 Nd 



(23.8%), 145 Nd (8.3%), 146 Nd (17.2%), 148 Nd (5.7%), 
and 150 Nd (5.6%). The lightest isotope is only produced 
in the s-process, the two heaviest isotopes are almost ex- 
clusively produced in the r-process, and the other four 
isotopes are produced in a combination of the two pro- 
cesses. The two odd- . A isoto pes exhibit HFS structure, 
like Sm. iRosner et al.l (2005) have measured the isotope 
shifts and HFS A and B constants for Nd transitions be- 
tween 4175 and 4650A. For the lines examined in our 
study, the Nd 11 HFS is comparable in scale to the Sm II 
HFS, < 0.05A, while the Nd 11 isotope shifts are rela- 
tively small, usually only < 0.025A. iDen Hartog et ail 
(2003) measured log(gf) values for many of these tran- 
sitions and noted that several of the yellow-red lines 
(mostly redward of 5100A) exhibited a doublet or triplet 
structure in their FTS profiles. 

We collect the Nd isotopes into fractions that might 
provide a good indication of pure-s- and pure-r- 
nucleosynthesis while making sensible choices about 
which sets of isotopes are shifted the most from one an- 
other. We adopt the Nd isotopic fraction defined by 

A^( 142 Nd) + iV( 144 Nd) 

/l42+144 iV(Nd) ' {} 

Predictions for this isotopic fraction are displayed in Ta- 
ble [H Note that the pure-s and pure-r predictions of 
the stellar model and classical method generally agree 
well for /142+144. In Appendix [Cl we present our calcula- 
tions of the HFS structure patterns for 6 Nd lines in the 
wavelength region covered by our spectra. 

4. ANALYSIS 

4.1. Radial Velocities 

To measure the radial velocities, we find the cen- 
troid of ~ 17— 19 Fe 1 lines and ca librate to the wave- 
lengths in ILearner & Thornd £[988) 10 , which are accu- 
rate to better than 0.3 mA. We find our wavelength 
scale zeropoint from a small number of telluric lines 
in the reddest orders of our spectra. Heliocentric ve- 
locities for each observation are calculated using the 
IRAF rvcorrect task. For HD 175305 and HD 122563, 
we measure heliocentric radial velocities of —184.2 ± 
0.3km s~ 4 and —26.0 ± 0.3km s" 1 , respectively, which 
are both consistent with measurements made by pre- 
vious studies to suggest no radia l velocity varia- 
tion (HD 175305: Car ney et all 2003; Nordstron Tet all 



2001 HP 122563: Bond 1985 IGratton fc Snedenl 1991 
Bar buv et alj 2003: Ho nda et al. 2004aHAokiet aljbood. 
20071) . For HD 196944, we measure a heliocentric radial 
velocity of —166.4 ± 0.3km s _1 . The scatter in our radial 
velocity measurements of this star from individual obser- 
vations over the very small date range 2006 June 12-18 
is consistent with a single value. iLucatello et al.l (2005) 
compiled recent radial velocity measureme nts for this 
star (ranging from —174.76 ± 0. 36 km s" 1 . | Aoki et al.l 
12001 to -168.49 ± 0.11km s" 1 , ILucatello et all 120051 . 



10 IWhaling et aD pOtjl ) propose a change in the 
ILearner~fc*?horng 1 119881 ) Fe I wavenumber scale, w ho cali- 
brated to the absolute Ar II wavenumber scale from INorlen I 
(1973). The proposed change, a multiplicative correction of 
1+67(8) xlO -9 , amounts to a correction of ~ 0.001 - 0.002 cm- 1 
or 0.3 — 0.5 mA in the visible spectral region examined in our 
study. This correction is negligible for our purposes. 
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see additional references therein) and concluded that it 
does exhibit radial velocity variations, consistent with 
their assertion that all carbon-enhanced, s-process-rich, 
very metal-poor stars are members of binary (or multi- 
ple) star systems. Our measurement supports the inclu- 
sion of HD 196944 in this category of stars. 

4.2. Equivalent Width Measurements 

Our high-resolution and high-S /N observations of these 
stars are superior to spectra analyzed in previous stud- 
ies, so we re-derive the atmospheric parameters for these 
stars. W e combine t he lin el ists and adop t log(gf ) val- 
ues from iFulbrightJ (|2000t >. IHonda et all (|2004al ). and 
llvans et al.l (|200a L Equivalent widths of isolated lines 
are measured in IRAF by fitting Gaussian profiles or 
by direct integration. We later discard the equivalent 
width measurements of lines that are approaching sat- 
uration, as indicated by their position on an empirical 
curve-of-growth. In each of our stars, equivalent widths 
are measured for ~ 40—50 Fe I lines, ~ 8 Fe II lines, 
~ 10 Ti I lines, and ^7 Ti n lines. A comparison of 
our equivalent width measurements with previous high- 
resolution studies of these reveals that there are no sys- 
tematic differences in the measured equivalent widths 
between our study and their studies. We find mean 
offsets of only A EW = (previous study) — (our study) = 
+0.5 ± 1.2 mA for HD 122563 (jHonda et all l2004aj ), 
AEW= -1.0 ± 2.2 mA for HD 175305 (|Fulbrightl[2000l) 
and A EW = -0.8 ± 2.3 mA for HD 196944 (jZacs et al.1 
fl99Sft . 

4.3. Atmospheric Parameters 

We use lKuruci ()1993f ) model atmospheres without con- 
vective overshooting for HD 175305 and HD 196944 and 
with overshooting for HD 122563 because its metallic- 
ity is beyond the range of atmospheres provided without 
overshooting. Interpolation software for the Kurucz grids 
has been kindly provided by A. Mc William and I. Ivans 
(2003, private communication). 

We have obtained V — J, V — H, and V — K col- 
ors for our stars from the SIMBAD 11 and 2MASS 12 13 
(jSkrutskie et al.1120061 ) databases. We adopt zero redden- 
ing values for these s tars as recommended by previous 
studies (HD 122563: IHonda et al.l I2004bb HD 175305 : 
Wallerstein et al.l 119791 : HD 196944: IZacs et all fl998l 



Jonsell et al.1 l2005f T We derive effective temperatures 



for these stars using the empiri c al col or-T e ff calibrations 
given by iRamfrez fe Melendezl (|2005f ). using interpola- 
tion software kindly provided by I. Ramirez 14 . Their 
work recalibrated and extended the range of color and 
metallicity applicability of earlier work by lAlonso et al.l 

11 This research has made use of the SIMBAD database, oper- 
ated at CDS, Strasbourg, France. 

12 This research has made use of the NASA/ IPAC Infrared 
Science Archive, which is operated by the Jet Propulsion Labora- 
tory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 

13 This publication makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the Univer- 
sity of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 



(|1996t Il999l ), who employed the infrared flux method to 
determine effective temperatures of F, G, and K dwarfs 
and giants. There is a rather small amount of scatter 
among giants in the temperature and metall icity ranges 
we are using (see IRamfrez fe Melendez|[2005T s Figure 4). 
The uncertainties intrinsic to the scatter of these calibra- 
tions and in the scatter of the temperatures predicted by 
each individual color lead us to adopt ±100K as the 
uncertainty in our effective temperatures. 

To determine the remainder of the model atmosphere 
parameters we use the most rec ent version of the LTE 
spectrum analysis code MOOG (jSnedenl [l973l ). We al- 
low for slight (~20K) modifications to the photometric 
temperatures to ensure that the abundances derived from 
lines with both high and low excitation potentials agree 
in our final atmospheric model. We derive the surface 
gravity by requiring that the abundances derived from 
neutral and ionized lines of Fe and Ti agree; we adopt 
±0.3 as the uncertainty in log(g). The microturbulence 
is measured by requiring that the abundances derived 
from strong and weak lines of each atomic species agree; 
we adopt ±0.3 km s _1 as the uncertainty in our micro- 
turbulence. These parameters arc varied iteratively until 
we arrive at our final values. 

We employ synthetic spectra to measure the macrotur- 
bulcnt broadening in our stars by fitting the line profiles 
of ~5— 10 clean Fe I and Fe II lines in each star. Un- 
certainties of 0.3— 0.4km s _1 were found for the macro- 
turbulent velocities. Our line profiles are well-fit by the 
convolution of the Gaussian (instrumental) and macro- 
turbulent (stellar) broadening terms. 

Our metallicity is defined to be the derived Fe I abun- 
dance, with an adopted uncertainty of ±0.15 dex. This 
uncertainty represents the systematic uncertainties as- 
sociated with determining T e g (~ 0.09-0.12 dex) and 
the intrinsic line-to-line scatter (~ 0.10 dex). Each of 
these stars are known to be q -enhanced (HD 175305: 
[a/Fe] = ±0.3. iFulbrightJ l25O0t HD 196944: [a/Fe] = 
+0.35. IZacs et al.l Il998t HP 122563: [a/Fe] = +0.4, 
IHonda et al.ll2004blh We have increased the overall metal 
abundance in our model atmosphere, [M/H], by 0.2 dex 
relative to the Fe I abundance to accoun t for the ex- 
tra e l ectron-donating q-e lements (see, e.g. iBrown et al.l 
119911: iSneden et~al"1 1 19941 ) 15 . For elemental abundance 
ratios (e.g., [X/Fe]) measured in our study, we reference 
neutral species to the Fe I abundance and ionized species 
to the Fe n abundance. 

Our final model atmosphere parameters are listed in 
Table [TJ For comparison, we show the atmospheres de- 
rived in previous studies of these stars. Table [T] indicates 
that our atmospheric parameters are in good agreement 
with these studies. 

4.4. Line Profile Analysis 

15 We have performed a test to compare abundances measured 
from equivalent widths of Ti, Fe, and Eu using three different Ku- 
rucz model atmospheres. For a typical metal-poor star, we adopt 
T cS = 5000 K, log(g) = 2.0, and v t = 2.0 km s -1 in all three mod- 
els. For atmosphere 1, [Fe/H] = -2.0 and [a/Fe] = +0.0. For 
atmosphere 2, [Fe/H] = —1.8 and [a/Fe] = +0.0. For atmosphere 
3, [Fe/H] = -2.0 and [a/Fe] = +0.4. Atmospheres 1 and 2 use the 
Kurucz grid of Q-normal atmospheres, while atmosphere 3 uses the 
Kurucz grid of a-enhanced atmospheres. The derived abundances 



are equivalent to within 0.02 dex in all cases, and this systematic 
14 Available online: https : //webspace .utexas . edu/ir68/tef f /dataliffctih is completely negligible for purposes of our study. 
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To measure the isotopic fractions, the observed absorp- 
tion lines are fit by synthetic spectra. Although our line 
lists are initial l y gen erated from the extensive lists of 
iKurucz fe Belli (|1995T l 16 . we employ experimental wave- 
lengths and \og(gf) values for all lines of interest and 
blended features whenever possible. Our Eu, Sm, and 
Nd lines are synthesized by accounting for the individual 
hyperhne and isotopic components; the isotope fraction 
and elemental abundance are varied to provide a best ht 
to the observed profiles. 

At the spectral resolution employed in this study, some 
of our unblended Fe I and II lines exhibit a slight pro- 
file asymmetry. This observation is consistent with our 
assumption that the macroturbulent velocity arises, at 
least in part, from large-scale convective motions in 
the stellar photosphere. These asymmetries account for 
w 1—2% of the total equivalent width of the line, and 
the typical bisector amplitude is « 200— 300ms" 1 . The 
shapes of our line bisectors are consiste nt with those ob - 
served for K giants, e.g., Arcturus, by iDravinU (|1987f ). 
We have not explored any attempts to introduce 3D 
modeling of the stellar atmospheres into our analysis; 
such models are not yet available for metal-poor giant 
stars. In this study we restrict ourselves to the use of 
one-dimensional, plane-parallel model atmospheres com- 
puted for the LTE case. 

In contrast to the challenge of detecting the presence of 
6 Li from the Li I 6707A line profile, which likely requires 
3D-NLTE modeling to account for the con vective asym- 
metry that can mimic the presence of 6 Li (|Cavrel et al.l 
2007), the change in the line profile that we analyze is 
very different. The presence of the 6 Li component may 
account for only a few tenths of a percent of the total 
flux in the red wing of the 7 Li line, which is a very subtle 
effect. The magnitude of the HFS for Eu causes gross 
changes in the shape of the entire line profile (see § 16. lj) . 
Sm and Nd each have seven naturally-occuring isotopes, 
so the shape of the line profile will be asymmetric along 
all parts of the line and not just along the red wing where 
convective asymmetries are most prominent. Therefore 
our choice to compute ID synthetic line profiles in LTE 
will not significantly affect our results for these species. 

4.4.1. Detailed Examination of the Absorption Line Profile 

Isotope shifts and HFS alter the shape of the line pro- 
file, and the isotopic fraction can be assessed by exam- 
ining the line profile from detailed comparisons between 
the observed and synthetic spectra. This method has 
be en used previously t o asse ss the Eu i s otopic fraction 
by ISneden et all (|2002l) and lAoki et all (|2003al lbl). We 
quantify our measurements of the isotopic fractions using 
a x 2 a l gorithm f after e.g..|Smith et al1l2001t lAoki et all 



l2003aHAsplund et al.ll200 



X 



(Oi - a 



(4) 



Here O — C is the difference between the observed and 
synthetic spectra at the ith spectrum point and a is de- 
fined as <Ji = [(S/N)x(/ ( ) 1/2 ]~\ where S/N is the ap- 
proximate signal-to-noise ratio of the continuum and / 
is the depth of the ith point relative to the continuum. 



Our best estimate of the isotopic fraction and abundance 
for each line is given by minimizing x 2 ■ Uncertainties 
are estimated at 1, 2, and 3cr confidence intervals from 
^X 2 — X 2 — Xmi7i — L 4, and 9, respectively. We con- 
sider the number of degrees of freedom in the fit to be the 
number of points in the spectral window plus the contin- 
uum normalization, abundance, and wavelength offset, 
all of which were allowed to vary in our fits. Reduced 
X 2 values are ~ 1, differing by less than a factor of two, 
which is still larger than would be expected for a nor- 
mal distribution of errors (e.g.. iPress et all H992). We 
attribute the scatter in the reduced x 2 values to our es- 
timation of the S/N ratio. The relative x 2 values for a 
given line, represented by A% 2 , however, do exhibit well- 
determined minima. We fix the macroturbulence, as de- 
termined in § 14.31 but we also estimate the uncertainties 
in the isotopic fraction introduced by the uncertainty in 
the macroturbulence, Vmacro ± AVmacro- Only our Eu 
and Sm lines with complete sets of hyperfinc constants 
(see § [3]) can be assessed by this method. 

4.4.2. Set the Absolute Wavelength Scale from Nearby 
Wavelength Standards 

In principle, the absolute wavelengths of the observed 
and synthetic spectra can by matched using wavelength 
standard Fe I and Ti n lines, and then the peak wave- 
length of the rare-earth line of interest can be measured. 
The isotopic fraction can be assessed because the iso- 
tope shifts — much more t han the HFS — af f ect th e peak 
absorption wavelength. iLundqvist et al.l (|2007f ) have 
used this method to qualitatively demonstrate that the 
Sm isotopic mix in CS 31082-001 is consistent with an 
r-process origin. 

For our Sm lines, a small wavelength shift 
(~ 0.02— 0.04A) is the dominant effect of changing the 
isotopic fraction; for our Nd lines, the shift is even smaller 
(~ O.OlA). In both cases, broadening of the line profile 
is a secondary effect. Therefore, given the small mag- 
nitudes of these shifts, it is critical that we match the 
wavelength scales of our synthetic and observed spec- 
tra with great accuracy. The absolute wavelengths of 
some of our Sm HFS components are accurate to ± 1 mA, 
as determined by the FTS measurements of 15 Sm by 
ILundqvist et al.1 120071 (see Appendix [B| , while the ab- 
solute wavelengths of the Sm lines not included in that 
study are accurate to « 2 mA. Unlike Sm, however, the 
atomic data for Nd does not include FTS measurements 
of individual isotope positions, so our reported absolute 
wavelengths of Nd components are based solely on mea- 
surements of the Nd energy levels. To avoid mixing 
measurements of energy level s from different stu dies, we 
adopt the energy levels from iBlaise et all (|1984l ). They 
estimate an uncertainty of ±0.005 cm" 1 for each level. 
Doubling this uncertainty translates to an uncertainty 
of - 1.5-2.0 mA across the region where our Nd lines 
are located. We note that the relative HFS component 
wavelengths — measured from FTS spectra — are more ac- 
curate than this in all cases. The abso lute wavelengths of 
the w avel ength standard li nes (Fe I: iLearner & Thorne 
1988 and Nav e et al.lll99l Ti II: iPickering et al.ll2001l 



20021 ) have a reported absolute accuracy in their FTS 
measurements of ~ ±0.001 cm" 1 (< 0.3 mA) in our 
wavelength regions. The uncertainties from our wave- 
16 Available online: http://kurucz.harvard.edu/iineiists.htmi length scale and co-addition of individual orders are no 
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larger than ± 0.9 mA, and we have verified that, after the 
wavelength solution has been applied, the Th-Ar emis- 
sion lines are correctly identified to within ~ ± 1 mA. 
All of these uncertainties are sufficiently small to enable 
reliable measurement of the isotopic fraction. 

In each echelle order containing a rare-earth line of 
interest, we locate 1-6 wavelength standard lines and 
match the wavelength scales of the observed and syn- 
thetic spectra at each of these lines, typically with a pre- 
cision of ~ 1 mA per line. When multiple wavelength 
standards are used, the scatter in our measured wave- 
length offsets necessary to ensure the matching of each 
line is ~ 5 mA, much greater than anticipated. This scat- 
ter appears to be independent of the ionization state, 
species, excitation potential, and relative strength of the 
wavelength standards, each of which could be expected 
to introduce small systematic offsets am ong lines with 
these different characteristics (s ee, e.g., iDravins et al.l 
119811 119861 : lAsplund et "all I2000L The isotopic fraction 
we would derive from our rare-earth lines of interest is 
therefore entirely dependent upon our choice of which 
wavelength standard line(s) we adopt and how many of 
these lines are present in the echelle order with the rare- 
earth line. 

To further understand this matter, we have also at- 
tempted to match the observed and synthetic spectra for 
wavelength st andards in th e Sola r spectrum. We com- 
pare both the iKurucz et al.1 (|1984h integrated-disk Solar 
spectrum, obtained with the FTS on the McMath So- 
lar Telescope at Kitt Peak National Observatory, and 
daylight sky spectra, obtained using the same instrumen- 
tal configuration on the McDonald Observatory 2dCoude 
echelle spectrograph as our stellar observations and re- 
duced in an identical fashion. It is well-known that astro- 
nomical echelle spectrographs and cameras can introduce 
a variety of non-linear distortions to the wavelength scale, 
whereas the spectrum obtained with the FTS should be 
immune to many of these distortions. We synthesize 
a spectral region from 4997-5010A (one echell e order) 
using an interpolated Hol weger fe M ullcr (1974) empir- 
ical Solar model photosphere. We apply our matching 
technique to the five unblended Fe I wavelength stan- 
dard lines in common to the two spectra. We measure 
the same line-to-line wavelength offsets for these lines in 
both spectra (agreement within our measurement preci- 
sion, ~ 1 mA) , suggesting that the wavelength scale of 
our 2dCoude spectra is consistent with the Kitt Peak 
FTS spectra. Nevertheless, the scatter in the absolute 
wavelength scale as determined from each of these lines 
in the FTS spectrum is still ~ 5 mA, indicating that there 
is a more fundamental problem that does not result from 
our spectral reduction techniques or our ability to match 
the observed and synthetic spectra. Unable to further 
identify and quantify the source of this discrepancy, we 
are unable to proceed with confidence by this method of 
analysis. 

4.4.3. Set the Relative Wavelength Scale at the Point of 
Insensitivity to the Isotopic Mix 

We also develop one additional method of quantify- 
ing the isotopic fraction that does not rely on nearby 
wavelength standard lines or fits to the shape of the line 
profile. The reason for developing such a method is our 
desire to assess the isotopic mix from all lines with max- 



imum isotopic shifts > 15 mA and from which an ele- 
mental abundance can be measured. This allows us to 
examine additional Sm lines that do not have full sets of 
measured hyperfine A and B constants as well as all of 
the Nd lines. 

This method attempts to match the observed and syn- 
thetic spectrum using the point in the synthetic spectrum 
that is insensitive to the s- or r-process mix. First, for 
each line, an isotopic fraction is assumed in the synthesis, 
and we always adopt / = 0.50 by default. Then an ele- 
mental abundance is derived from a synthesis with this 
isotopic fraction. Using this elemental abundance, two 
new syntheses (one with a pure-s-process mix of isotopes 
and one with a pure-r-process mix) are created. The 
intersection of these two syntheses is defined to be the 
point of insensitivity to the isotopic fraction. The relative 
wavelength offset between the point of insensitivity and 
the observed spectrum is then adjusted until the point 
of insensitivity is coincident with the flux of the nearest 
pixel of the observed spectrum. Next, we adopt this rel- 
ative wavelength offset and generate a set of synthetic 
spectra with a range of isotopic fractions but the same 
elemental abundance. A new isotopic fraction is then 
measured by minimizing the sum of the squares of the 
differences between the observed and synthetic spectra. 
This isotopic fraction is then used as the input to derive 
an elemental abundance, and the process is repeated un- 
til the elemental abundance and isotopic fraction have 
each converged. We show an example where the point of 
insensitivity has been matched to the observed spectrum 
in Figure [3 

The uncertainty in the matching process is given by 
the width of one pixel; therefore the uncertainty in 
the derived isotopic fraction is estimated by shifting 
the synthetic spectrum relative to the observed spec- 
trum by ±U!pj xe i/2. Even at a spectral resolution of 
R ~ 120, 000, the width of one pixel is still a substantial 
fraction of the width of the absorption line itself. This 
method also places a disproportionate amount of weight 
on the one or two flux measurements closest to the point 
of insensitivity, rather than taking account of the infor- 
mation from the full width of the line to match the ob- 
served and synthetic spectra. The precision achieved by 
this method is hardly adequate to clearly assess the iso- 
topic mix — particularly for Nd — but it does provide a 
consistent set of measurements for all of our lines of Sm 
and Nd. 

5. NEUTRON-CAPTURE ELEMENTAL ABUNDANCES 

Our mean elemental abundance of Eu in HD 175305, 
[Eu/Fel = +0. 46 ± 0.07, is in good agreement with 
lFulbrightl(l2000l)' s abun dance, [Eu/Fe] = +0.48 ± 0.105, 
and lCowan et al. (2005)'s abundance, [Eu/Fe] = +0.55 ± 
0.10. Our synthesis of the Eu n 4205A line in HD 122563 
does permit measurement of the Eu elemental abun- 
dance, [Eu/Fe] = -0.55 ± 0.06, assuming 12 C / 13 C= 5 
(|Lambert fe Snededfl977t iHonda et ai1l2004bh and a Eu 
isotopic fraction of /151 = 0.5. The Eu abundance 
changes by ±0.03 dex when the assumed isotopic ratio 
is varied by ± 0.2 dex. This elemen tal abundanc e agree s 
very well with the value found by IHonda et al.l |2006) , 
[Eu/Fe] = —0.52, and with many previous studies of this 
star. The 4205A line in HD 196944 is far too sensitive 
to the C abundance to derive a meaningful Eu elemental 
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abundance. 

We measure Sm elemental abundances of [Sm/Fc] = 
-0.56 ± 0.12 from 1 line in HD 122563, [Sm/Fe] = 
+0.37 ± 0.11 from 9 lines in HD 175305, and [Sm/Fe] = 
+0.88 ± 0.08 from 6 lines in HD 196944. Previous stud- 
ies of Sm elemental abundances in these sta rs have found 
[Sm/ Fe] = -0.40 ± 0.17 in HD 122563 (|Honda et all 
[2001 and rSm/ Fel = +0.78 ± 0.23 in HD 196944 
(jAoki et"al"1l2002D ; Sm has not been measured previously 
in HD 175305. Unlike the uncertainties in the isotopic 
fractions, the uncertainties in the elemental abundances 
are dominated by the systematic uncertainties in the 
model atmosphere parameters and are affected somewhat 
less by internal scatter from the log(gf) values. 

We measure Nd elemental abundances of [Nd/Fe] = 
+0.22 ± 0.13 from 12 lines in HD 175305 and [Nd/Fe] = 
+0.95 ± 0.13 from 11 lines in HD 196944. Previ- 
ous studies of Nd elemental abu ndances in these s tars 
have found [Nd/Fe] = +0.34 (|Burris etal.1 l200l in 
HP 1 75305 and [Nd/Fe] = +0.93 ± 0.09 (iVan Eck et all 
l200lh . [Nd/Fe] = +0.94 ± 0.17 dZacs et al.l ll998D. and 
[Nd/Fc] = +0.86 ± 0.20 (|Aoki et aTl 2002) in HD 196944. 

We have also measured elemental abundances for five 
additional heavy n-capture elements in HD 175305. 
These abundances are listed in Tabic [3] along with our 
Eu, Sm, and Nd abundances. The uncertainties listed in 
the table reflect the uncertainties in the atmospheric pa- 
rameters, continuum placement, line-to-line scatter, and 
log(<7/) values, all added in quadrature. In cases where 
only one line of a species was available for analysis, we 
substitute the uncertainty in the spectral synthesis fit to 
the observed spectrum in place of the line-to-line statis- 
tical scatter in this formula. In Figure [3] we plot our el- 
emental abunda n ces, combined with m easu rements from 
Fulbri£h3 (|2000D . ICowan et alj (|2005h . and lLawler et all 
<l2007ft . ag ainst the scaled S. S . r-process abundance 
pattern of ISimmerer et all (|2004h . While not a perfect 
fit, this figure suggests that the scaled S. S. r-process 
abundance pattern provides a respectable fit of the ob- 
served n-capture abundances beyond Ba. HD 175305 
may not be a star with pristine r-process material, but 
the r-process appears to have contributed the majority 
of the n-capture species present. We have also used the 
Keck I HIRES spec trum of HD 175305, obtained for the 
ICowan et al.l (|2005l ) study, to determine an upper limit 
for Pb from the 4057A line. 

6. NEUTRON-CAPTURE ISOTOPIC FRACTIONS 

We emphasize that measurements of isotopic fractions 
are highly insensitive to the choice of atmospheric param- 
eters. Separate isotopic components of a line all share 
the same initial and final energy levels (since the HFS 
of these levels is negligible relative to the energy differ- 
ences themselves), eliminating systematic errors associ- 
ated with the excitation and ionization states of each 
isotope. 

6.1. The Europium Isotopic Fraction 
6.1.1. In the Sun 

For S. S. material, the n-capture isotopic abundances 
are precisely known from studies of terrestrial meteorites. 
The Eu II i sotopic f ractio n in t he Solar spect rum was first 
studied by iHaugd (|1970l ) and|HauggJ (| 19721 ). who found 



i53 Eu/ i5i E = (48 ± 6)/ ( 5 2 + 6) (or f 151 = 0.52 in our 
notation). lLawler et al.l (|2001bl ) reevaluated the Solar 
Eu isotopic fraction, finding /151 = 0.50 ± 0.07 from the 
4129A line. We do not repeat this measurement here. 

6.1.2. In Metal-Poor Stars 

From the Eu 11 4129A line in HD 175305 we measure 
an isotopic fraction of /151 = 0.50 ± 0.04, as shown in 
Figure [U The A x 2 values are shown in Figure O for dif- 
ferent values of the isotopic fraction and macroturbulent 
velocity. The minimum A \ 2 values have well-defined 
minima, leading to the small uncertainties in the derived 
isotopic fraction. As noted in Appendix lA.il we cannot 
trust our syntheses of the Eu 11 lines at 4205 and 4435A 
at the precision necessary to measure an isotopic frac- 
tion. For the sake of comparison, in Figures [5] and [7] we 
show our synthesis of these two lines in HD 175305 given 
our best estimate of the blending features, but we em- 
phasize that these syntheses display the best-fit isotopic 
fraction derived from the 4129A line. 

The Eu 11 4129A line was not covered in our spectra of 
either HD 196944 or HD 122563. 

Although only two naturally-occurring isotopes of Eu 
exist, only a small spread is predicted between f( 51 
and /C 51 (see panel a of Figured]). It is difficult even 
with very high-resolution (R > 100, 000) and high- 
S/N (~ few hundred) data on one or more lines to 
achieve sufficiently small uncertainties on isotopic abun- 
dances to determine t he nucleosynthetic source of the Eu. 
ISneden et all (|2002f ); lAoki et all (|2003aU bT) have success- 
fully shown that the Eu s- and r-process nucleosynthetic 
signatures can be distinguished at the isotopic level, us- 
ing approximately three lines in each of seven stars. The 
measurements made from individual lines in each star in 
these studies have a large degree of scatter, causing the 
final Eu isotopic fraction to be mildly sensitiv e to the set 
of lines chosen. For example, we note that lAoki et all 
l|2003al) 's two measurements of /151 from the 4435A line 
appear systematically higher than the results obtained 
from the 4205A line by « 10-15%, which is roughly the 
same as the change in the Eu isotopic fractions that are 
trying to be distinguished. Although we do not measure 
an isotopic fraction from the 4435A line in HD 175305, 
our syntheses hint at a similar result in the small region 
without obvious blends (4435. 40-4435. 50A). Additional 
measurements of the Eu isotopic fraction from these lines 
would be helpful to confirm or refute these possible sys- 
tematic trends. If such systematic trends exist, they 
likely result from the presence of blending features, and 
the measured isotopic fractions may be as sensitive to 
the blending features as to the Eu isotopes themselves. 
There is little hope of being able to distinguish Eu iso- 
topic fractions with any greater precision than has al- 
ready been achieved using these methods on this set of 
lines, including the Eu 11 lines at 3819 and 3907A. 

6.2. The Samarium Isotopic Fraction 
6.2.1. In the Sun 

lEkeland fc Haugel (|1975f ) analyzed the Sm isotopic 
fraction in the Solar photosphere from two lines, 4467 
and 4519A, finding isotopic fractions consistent with the 
meteoritic values. The 446 7 A line is broadened slightly 
by HFS, with almost zero broadening resulting from the 
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isotope shifts, while the 4519A line is broadened slightly 
by isotope shifts. High-precision single- frequency laser 
measurements of the isotope shifts for 147 Sm and 149 Sm 
have not been measured in the lab for the 4519A line, 
therefore we do not attempt synthesis of this line. 

We synthesize the 446 7 A line using the 
iHolweger fe Mullerl (I1974T ) model atmosphere and 
the iKurucz et alj (|1984D Solar spectrum. We fix the 
broadening from nearby Fe lines. The total spread in 
the HFS and isotope shifts for this line are < 0.05A 
and < O.OlA, respectively, and we cannot measure the 
isotopic fraction, either in terms of / oc id or /152+154, 
from the line-profile-fitting m ethod even when a high - 
quality Solar spectrum is used. lEkeland fe Haugel (|1975f ) 
provide very little information about their method of 
measuring the Sm isotopic fraction, and their very 
small uncertainties (/odd = 0.28 ± 0.07) for a line with 
relatively small HFS and isotope shifts give us some 
reservation in accepting their result. All of the other Sm 
transitions listed in App endix I A . 2 1 ar e either far too weak 
or are severely-blended with strong transitions of other 
species to measure reliable isotopic fractions in the Solar 
spectrum, so we are unable to make any measurements 
of the Sm isotopic fraction in the Solar photosphere. We 
remind the reader that the blended lines may be reliable 
for measurements of the Sm elemental abundance, but 
we cannot trust our syntheses of the blending features 
at a precision necessary to measure isotopic fractions. 

6.2.2. In Metal-Poor Stars 

Figures [5] and [5] show our syntheses of the Sm 11 line 
at 4424A in HD 175305 and HD 196944, respectively, 
and the results from our line profile fits. We also in- 
clude the Cr abundance of the Cr 1 4424. 28A line, which 
blends the blue wing of the Sm line, as a free parame- 
ter. In HD 122563, this Sm line has a relative depth of 
~2%. We do not claim measurements of either the el- 
emental abundance or isotopic fraction of Sm from this 
line in HD 122563; it is clear that this line is simply too 
weak to warrant further analysis. In HD 175305, we mea- 
sure a Sm isotopic fraction of /152+154 = 0.55 ± 0.14. 
In HD 196944, we measure a Sm isotopic fraction of 
/152+154 = 0.35 ^0 20- In addition to showing the best fit 
curves and their uncertainties in Figures [8] and [H we also 
show the syntheses for pure-s- and pure-r-process nucle- 
osynthesis and the (O — C) curves for these syntheses, 
illustrating the contrast between the two extremes. In 
each of these figures, the pure-s- and pure-r-process syn- 
theses in panel b are compared to the observed spectra 
in panel a, illustrating the degree to which the observed 
spectrum tends toward one synthesis or the other. 

FigureflOlshows our synthesis of the Sm 11 line at 4604A 
in HD 175305. We measure a Sm isotopic fraction of 
/152+154 = 0.501q °4. This line is weak, with a relative 
depth of only 5%. Figure QT] shows the A% 2 values for 
our fit to this line; A \ 2 values from the Sm 4424A line 
in both stars are similar to the 4604A line. 

If instead we match the observed and synthetic spec- 
tra at the point of insensitivity to the isotopic frac- 
tion, we find /i52+i54 = 0.41 ± 0.32 from the 4424A 
line, /152+154 = 0.49 ± 0.17 from the 4604A line, 
and /152+154 = 0.49 ± 0.20 from the 4719A line in 
HD 175305. In HD 196944, we find /152+154 = 0.37 ± 



0.28 from the 4424A line. While we regard the line profile 
fits as the superior measurements of the isotopic fraction, 
the agreement between the sets of measurements from 
the two methods is reassuring. 

6.3. The Neodymium Isotopic Fraction 

6.3.1. In the Sun 

We check our Nd isotopic fraction measurements 
against the Solar spectrum. We are aware of no pre- 
vious attempts to examine the Nd isotopic mix in the 
Solar photosphere, but the individual S. S. isotopic abun- 
dances have been established from studies of terrestrial 
meteorites. The Nd 11 line at 4446A would appear to 
be the most reliable indicator of the Nd isotopic fraction 
due to its relative strength and freedom from blending 
features. The line profile shape does not change enough 
to allow us to apply our fitting algorithm to measure 
the isotopic fraction, and we can find no good candidate 
Fe 1 or Ti 11 lines nearby to match the observed and syn- 
thetic spectra. By using the method of matching the 
observed and synthetic spectra at the point of insensi- 
tivity to the isotopic mix, we derive an isotopic fraction 
°f /142+144 = 0.53lg53, which is unfortunately unable 
to exclude any possible combination of isotopes. There- 
fore, we are unable to provide independent confirmation 
of the S. S. Nd isotopic fraction from analysis of the Solar 
photosphere. 

6.3.2. In Metal-Poor Stars 

Because the isotope shifts of Nd are smaller than Sm 
and because we cannot match our observed and syn- 
thetic spectra with a high absolute accuracy, we can 
only rely on the third matching method to measure 
the isotopic fraction of the Nd lines in these stars. In 
HD 175305, we measure /142+144 = 0.25^25 from the 
4177A line, /142+144 = 0.1 from the 4446A line, 

and /142+144 = 0.33tjj;g from the 4567A line. The 
mean isotopic fraction of these three lines is /142+144 — 
0.2llg 2p Combining these measurements together us- 
ing the asymmetric uncertainties that reflect the fact that 
the isotopic fraction cannot be less than 0.0 or greater 
than 1.0 could introduce a bias into the weighted mean. 
Therefore we have estimated the uncertainties in these 
measurements also by assuming (non-physically) that the 
isotopic fraction can exceed these limits. We find un- 
certainties of ±1.10, ±0.85, and ±1.00 on these three 
lines, respectively, resulting in a weighed uncertainty of 
±0.56 on the mean. Our final quoted uncertainties do 
impose the physical limits on the isotopic fractions. In 
HD 196944, we measure /142+144 = 0.36+^;^ from the 
4446A line. The non-physical uncertainty on this mea- 
surement is ± 0.68, which provides a relative comparison 
to our ability to measure the 4446 A line in HD 175305. 
Unfortunately none of these individual measurements 
provides any real insight into the Nd isotopic fraction 
due to the size of the uncertainties. The large uncertain- 
ties also prevent us from pointing out even a contrast 
between the Nd isotopic mixes in the two stars, much 
less make any detailed assessment of the s-process nucle- 
osynthetic predictions. 

7. INTERPRETATION OF THE NUCLEOSYNTHETIC 
SIGNATURES 
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In HD 175305, our Sm isotopic fraction /152+154 = 
0.51 ± 0.08, a weighted average of our measurements 
from profile fits to the 4424 and 4604A lines, suggests 
an r-process origin. It is tempting to surmise that our 
Nd isotopic fraction, /142+144 = 0.2lt°;^, derived fr om 
three lines, is suggestive of an r-process origin, but the 
large uncertainty cannot exclude an s-process origin. 
Our Eu isotopic fraction, /151 = 0.50 ± 0.04, derived 
from only the 4129A line and shown to be a reasonable 
match for the 4205 and 4435A lines , implies an r-process 
origin by the ICowan et al.l (|2006l ) predictions but ex- 
clu des neither a pu r e-s- n or a pure-r-process origin by 
the lArlandini et al.l (|1999l ) predictions. It is interesting 
to note that our isotopic fractions in HD 175305 suggest 
a nucleosynthetic history predominantly — but perhaps 
not purely — dominated by the r-process, which also is 
suggested by the elemental abundance trends. The fact 
that both the elemental and the isotopic distributions in 
HD 175305 agree with the predicted elemental and iso- 
topic distributions of r-process material in the S. S. adds 
to the preponderance of evidence that supports the hy- 
pothesis of a universal r-process mechanism for elements 
with Z > 56. 

In HD 196944, our measured Sm isotopic fraction, 
/152+154 = 0.35 ± 0.14, is suggestive of an s-process 
origin. Our Nd isotopic fraction, /142+144 = 0.36+0 36 
can exclude no possible values of the isotopic mix. The 
interpretation for the Sm is consistent with the s-process 
elemental abundance signatures and radial velocity vari- 
ation described in § 14.11 

Our best measurements for each species are summa- 
rized in Figures [12] and [13] for both HD 175305 and 
HD 196944. In Figure [HI the differences between the 
pure-r- and pure-s-process predictions are scaled to- 
gether to emphasize the relative precision with which 
the two processes can be distinguished by each species. 
In contrast, Figure [TBI displays our measurements in an 
absolute sense to emphasize the relative precision with 
which each species was measured. Expressed this way, 
it is clear that Eu can be measured precisely but does 
not easily reveal its nucleosynthetic origin because of the 
small difference between /f 51 and /f 51 . The uncertain- 
ties on each Nd line can hardly exclude any possible value 
of the isotopic fraction. Only the Sm isotopic fractions, 
which have uncertainties on the individual measurements 
that are smaller than the difference between /[52+154 and 
/i52+i54> on?er anv opportunity to distinguish the two 
processes with reasonable precision. 

8. EU/PB AND 151 EU/PB ABUNDANCE RATIOS 



In § 13.21 we pointed out the great disc r epanc y be- 
tween the stellar mode l of lArla ndini et al (1999|) and 
the classical method of ICowan et al.l ([2006 ) in predict- 
ing the amount of 151 Eu produced by the s-process: 
N s = 0.00304 ± 0.00013 and N s = 0.0001, respec- 
tively. The stellar model predicts that the s-process 
should produce ~ 30 times more 151 Eu than the classical 
method does. For a certain amount of 151 Eu produced by 
an s-process nucleosynthesis event, the classical method 
would then predict that the isotopes near the termination 
of the s-process — such as Pb — should be overproduced 
by a factor ~ 30 relative to the stellar model. 

We choose Pb for this comparative analysis because 
it is so dominantly produced by the s-process. In Pb- 



enriched metal-poor stars, such as HD 196944, the Pb 
overabundance often will reach [Pb/Fe] > +2.0, mak- 
ing it easy to detect at the Pb 1 4057A transition. Be- 
cause Eu is produced more easily by the r-process than 
by the s-process (relatively speaking) and Pb is produced 
more easily by the s-process than by the r-process, the 
Eu/Pb and 151 Eu/Pb abundance ratios in metal-poor 
stars should be extremely good discriminators of the pre- 
d icted abundance s of the se species. 

lArlandini et al.l (|1999f) do not predict pure-r-process 
Pb abundances due to the fact that Pb can be over- 
produced relative to the lighter-s-process nuclei in low 
metallicity AGB stars. In other words, at lower metal- 
licity a greater number of neutrons is captured per Fe- 
peak seed than at higher metallicity, so the s-process 
abundance distribution at low metallicity fav o rs th e 
heavest s-process nuclei ; see, e.g., iGallino et al.l (fl99l . 
IBusso et al.l (|1999h . and lStraniero et al.l (|2006f l 

In Figure Q3] we have compiled the derived Eu/Pb 
ratios (or upper limits) for 35 metal-poor stars from 
a variety of recent studies and for HD 175305. The 
Eu/Pb ratio shows no obvious trends when expressed 
as a function of [Fe/H], but we note that this sam- 
ple of Eu and Pb abundances is far from complete. 
Adopting the original authors' classifications for the 
enrichment patterns for these star s, this sample con- 
sists of 6 r-process-enriched stars (ICowan et al.l 1 2002 ; 
Hill et all [2001 lAoki efall l2003al: iSneden et all E)03; 
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A. Frebel, 200 7, private communication), 15 s- process- 
enriched stars (lAoki et al.ll2001l l2002t I Johnson fc Bolte 
| 2002t|Cohen et al.ll2003tlBarbuv et al.ll2005l : ICohen et al. 
l2006 HPreston et al.ll2006D and 5 (r + s) -pr o cess-enriched 
stars (llvans et al.1 120051: lAoki et al.l l2006t Uonsell et all 
l2006t iMasseron et al.ll2006[ ) from the Galactic halo pop- 
ulation. The r-process-enhanced sample also includes 
10 stars fro m the metal-poor globular clusters M 13 and 
NGC6752 (jYong et al.ll2006T) . The small number of r- 
process enriched stars results from the observational dif- 
ficulty of detecting Pb in stars without s-process en- 
richment. The Eu/Pb ratios of the r-process-enhanced 
halo stars and globular cluster stars all suggest that the 
r-process predictions are correct. Despite the fact that 
log £ (Eu/Pb) s differs by ~ 0.6 dex between the stellar 
model and classical method predictions, the measured 
abundance spread in the s-process-enriched stars is still 
too great to unambiguously distinguish the two predic- 
tions from this information alone. The s- and (r + s)- 
enhanced classes of stars show little difference from one 
another when their Eu and Pb abundances are expressed 
this way. 

In Figure [15] we have complied a subset of 7 of these 
stars with measured Eu isotopic fractions, and we plot 
log e ( 151 Eu/Pb) as a function of metallicity. This sub- 
set includes 5 r-process-enhanced stars and 2 (r + s)- 
process-enhanced stars. Again, the r-process-enhanced 
halo stars cluster around the r-process predictions. Here, 
though, the stellar model and classical method predic- 
tions for log e( 151 Eu/Pb) differ by ~ 1.7 dex, and the 
two r+s s tars with measured Eu isotopic fractions clearly 
favor the lArlandini et al.l () 19991 ) predicted 1 51 Eu abun- 
dance. This is not unexpected, given that A oki et al.l 
(|2003al) , who measured the Eu isotopic fraction in these 
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two stars, repo r ted t hat /151 was consistent with the 
lArlandini et al.1 (|1999h predictions. 

It is worth considering whether these r + s stars 
may have enough r-process enrichment to give them 
a log e ( 151 Eu/Pb) abundance that lies between the 
ICowan et all (poll pure- s and pure-r predictio ns, in- 
advertently suggesting the lArlandini et all (|1999f ) value. 
To investigate this matter, in Figure [16] we plot the 
log e (Ba/ Nd) abundance ra tios for the full sample of 
halo stars; lYong et"aH (|2006f ) did not measure Ba or Nd 
abundances in their globular cluster stars. While there is 
still a relatively large amount of scatter among the s- and 
(r+s)-enriched stars, neither population appears distinct 
from the other, and the two r+s stars shown in Figure [TBI 
scatter above and below the predicted pure-s-process 
abundance value. We therefore conclude that these 
two stars are representative of predominantly-s-process- 
enriched material, and these ob servations thereby sup- 
port the lArlandini et al.l ([1999) 151 Eu abundance pre- 
diction. We caution that this conclusion hinges on the 
Eu isotopic fraction measured in only two stars, and ad- 
ditional observations of the Eu isotopic fraction in s- 
process-enriched stars would be most helpful to further 
r esolve the matte r . 

iStraniero et al.l (|2006h show that in low-metallicity 
stars in the TP-AGB phase the abundance ratios of nu- 
clei near the first and second s-process peaks are only 
weakly affected by the efficiency of the 13 C pocket (and 
hence the neutron source for the s-process), whereas 
abundances of species near the third s-process peak are 
strongly dependent on this parameter. Therefor e the 
abundance predictions of the lArlandini et al. (1999]) stel- 
lar model, which was computed for ~ Solar metallicity 
stars in the TP-AGB phase, should predict the lighter 
s-process abundances at low metallicty better than it 
predicts the Pb abundance at low metallicity. In a low- 
metallicity s-process nucleosynthesis event, we would ex- 
pect the 151 Eu/Pb ratio to decrease relative to a higher- 
metallicity case as a higher fraction of seed nuclei are con- 
verted to Pb. Therefore, the stellar model pure-s Eu/Pb 
and 151 Eu/Pb ratios should represent upper limits for 
these values. In both cases, decreasing these ratios would 
bring the stellar model predictions into better agreement 
with the classical method predictions. 

9. PROSPECTS FOR FUTURE STUDIES OF 
NEUTRON-CAPTURE ISOTOPIC FRACTIONS 

Previous studies have demonstrated that the Eu and 
Ba isotopic fractions can be measured, we have shown 
that Sm isotopic fractions can also be quantitatively mea- 
sured, and we also suggest that Nd isotopic fractions can 
be measured, albeit with lower precision that the other 
rare-earth species. The measurements of each set of iso- 
topes give an indication of the stellar nucleosynthetic his- 
tory that is consistent with other indicators, and we now 
discuss the role that measurements of Ba, Nd, Sm, Eu, 
and other n-capture isotopic fractions might play in fu- 
ture studies. 

9.1. Isotopic Fractions of Multiple Neutron- Capture 

Species 

It is a success of both the stellar model and classical 
method that these models' predictions of the Sm /152+154 
isotopic fraction are in good agreement. Suppose that 



both Sm isotopic fractions /152+154 and / oc id could be 
measured from the same set of lines. This would en- 
able a comparison of the models' predicted abundances 
of 147 Sm, which differs by a factor of ~ 3 and affects 
the predicted / oc id values, whereas both models generally 
agree on 149 Sm, 152 Sm, and 154 Sm. One could measure 
the isotopic fractions of Sm and Ba in an r-enriched star. 
If a convincing argument could be made that the Ba and 
Sm should have originated from the same nucleosynthetic 
event(s) (i.e., an extra contribution has not been made 
to the Ba abundances by another n-capture process), the 
isotopic fractions measured for both Ba and Sm should 
show the same relative proportion of r-process isotopes. 
This scenario could permit a direct test to the predicted 
Ba isotopic fractions, which differ by 26% in their r-only 
predictions for /odd, with the most significant differences 
arising in the models' predictions for 137 Ba and 138 Ba. 
Such tests would be greatly welcomed. 

It is not obvious why Ba abundances are sometimes 
larger than the scaled S. S. r-process abundances (scaled 
to the heavier n-capture elements) in metal-poor stars 
with an r-only chemical signature and no evidence of bi- 
narity, such as HP 175305. (See al so, e.g., Figure 15 of 
iHonda et al.ll2004al or Figure 15 of iBarklem et al.ll2005l 
wher e some o f their r-enriched stars exhibit this behav- 
ior.) lAsplu nd (2004) suggests that the 3D Ba abundance 
correction may be as much as —0.5 dex for some metal- 
poor dwarfs and subgiants. While this matter may in- 
deed be a result of poor Ba abundance determinations 
due to ID LTE modeling or to uncertainties in the S. S. 
abundance predictions, the search for an astrophysical 
explanation for this phenomenon should be explored fur- 
ther. Measurements of Ba and other n-capture isotopic 
fractions could be used to identify any correlations that 
may exist between super-r-process Ba abundances and 
other nucleosynthetic signatures. 

The combination of Ba, Nd, Sm, and Eu isotopic frac- 
tions could provide more complete knowledge of abun- 
dances at the isotopic level over a range of ~ 20 mass 
numbers (A =134—154). This information would com- 
plement elemental abundance measurements of lighter 
and heavier n-capture species to get a more complete 
picture of the n-capture nucleosynthesis, constrain the 
conditions (e.g., temperature, neutron density, etc.) that 
are required to produce the r-process elements, and de- 
termine the actual r-process path by identifying the in- 
dividual isotopes that part icipate in this process (see 
ICowan &: Thielema nn 2004 for further discussion). Ex- 
tended discussion of the relationship between neutron 
densities in the "main" and "weak" r-processes and the 
isotopic abundance si g natur es in this mass range is pro- 
vided in iKratz et al.l (|2007l ) and will not be repeated 
here. There is evidence that different mass regimes may 
be synthesized by different astrophysical processes in ad- 
dition to the s- and r-processes, such as the light ele - 
ment primary process, or LEPP (Travag rio et al.ll2004h . 
While the LEPP would have the greatest effect on nucle- 
osynthesis of lighter elements such as the Sr-Y-Zr group, 
knowledge of isotopic data for elements with very differ- 
ent Z would still be extremely informative. 

This study has made use of the most current and com- 
plete laboratory atomic data to compute the isotope 
shifts and HFS patterns for Sm and Nd. While we are 
pleased that the available atomic data permits measure- 
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ment of the isotopic fractions of these two species, com- 
plete sets of HFS constants and isotope shifts have not 
been measured for many lines that lie in the yellow-red 
regions of the spectrum. These transitions are gener- 
ally weaker in stellar spectra, but the isotope shifts may 
be large enough nonetheless to permit reliable measure- 
ments of the isotopic fraction even with weak lines. We 
encourage such efforts from the atomic physics commu- 
nity in the future. We note that these lines may be 
significantly stronger in more metal-rich stars or those 
which exhibit enhanced abu ndances for s-process e le- 
ments, such as Ba stars fe.g.. lAllen fc Barbuvll2006li bl). 
where the blue lines we have used in our study are likely 
to be heavily blended with other atomic and molecular 
features; in these cases, lines in the yellow-red spectral 
regions may provide the only opportunity to establish 
the n-capture signatures at the isotopic level. 

9.2. Individual Advantages and Disadvantages of 
Europium, Samarium, and Neodymium Isotopic 
Fractions 

One advantage of using measurements of the Eu iso- 
topic fraction as a nucleosynthesis probe is the wide HFS 
structure that can be measured with only moderately- 
high spectral resolution (R ~ 45-60,000). In r-enriched 
stars, the Eu lines are strong enough to facilitate mea- 
surement of the isotopic fraction even at very low metal- 
licities. But, as we have argued in § 16.11 it is unlikely 
that Eu isotopes can be measured with the precision to 
unambiguously distinguish s- and r-process nucleosyn- 
thesis using the lines in the near-UV and blue regio ns of 
the spectrum if indeed the lArlandini et alj (|1999h pre- 
dictions are correct. 

Sm presents a different set of circumstances. The 
HFS of Sm is not nearly as wide as Eu, and the lines 
with the widest structure are also very weak (or un- 
detectable) and sometimes blended in very metal-poor 
stars. Though we are forced to parametrize the s- and 
r-process nucleosynthetic predictions into some set of iso- 
topes (e.g., / dd or /152+154) to facilitate distinguishing 
among seven isotopes simultaneously, the greater differ- 
ence predicted between /f 52 +i54 and /152+154 allows one 
to cleanly distinguish these different nucleosynthetic sig- 
natures with reasonable measurement uncertainties. In 
addition, if Sm isotope fractions could be measured in a 
large sample of stars, it may be possible to observe the 
changing contributions of s- and r-process nucleosynthe- 
sis in_star^_ovBr_rtie_liistOTy_^of the galaxy (cf. Figure 6 
of iMa shonkin a &: Zhaol l2006l who performed a similar 
study using Ba isotopes). Therefore, while Sm lines that 
permit measurement of isotopic fractions may not be ac- 
cessible in all metal-poor stars, we propose that in favor- 
able cases Sm isotopes may give a clearer picture of the 
n-capture nucleosynthetic history than Eu isotopes can. 

The Nd isotopic fraction is more difficult to measure re- 
liably than the Sm isotopic fraction. The smaller isotope 
shifts of Nd translate directly into larger uncertainties 
in the measured isotopic fraction on any individual line. 
This necessitates that as many clean lines as possible 
be measured in order to reduce the uncertainties on the 
mean Nd isotopic fraction. Only then is it possible to un- 
ambiguously distinguish distinct chemical signatures in 
the Nd isotopes. Nd is produced more by the s-process 
than the r-process, though, so Nd lines may be easier 



to observe in stars with an s-process chemical signature 
than either Eu or Sm. 

Of course these suggestions only apply to the examina- 
tion of the nucleosynthetic history at the isotopic level. 
At the elemental level Eu should continue to serve as an 
excellent barometer of r-process nucleosynthesis, and Ba 
or La should continue to serve as barometers of s-process 
nucleosynthesis. 

9.3. s-process Yields from Models of Low-Metallicity 

Stars on the AGB 

While the classical method of modeling the s-process 
is, by definition, model- independent, it has been shown 
to predict an overproduction of 142 Nd relative to other 
isotopes of Nd and pure-s nuclei, e.g., a s determined 
from studies of the Mu rchison meteorite (|Zinner et al.l 
ll99lUGuber et al.|[l997t) . This isotope is only produced 
by the s-process and contains 82 neutrons, one of the 
magic neutron numbers, making this nucleus relatively 
more stable than other nuclei with similar mass numbers 
and thus violating the classical method's assumption of a 
s moothly- varying o ~ N s cu rve. The complex stellar model 
of lArlandini et all (1999) can avoid this difficulty near 
the magic neutron numbers. 

Of the much theoretical and computational work that 
has been performed in recent years to better understand 
the structure and evolution of low-metallicity stars on 
the AGB, one goal has been to provide a complete set 
of chemical yields from H to Bi, at the termination of 
the s-process . The first such set of yields was presented 
bv ICristallol ((2006), for 2M Q stars on the AGB with 
Z = 1.5 x 10" 2 , 1 x 10" 3 , and 1 x 10" 4 . Similar mod- 
els surely will follow, and the observational verification 
for these models will rely heavily on isotopic abundances 
as well as elemental abundances, particularly for magic 
neutron number nuclei such as 142 Nd. The interpreta- 
tion of n-capture nucleosynthetic signatures of the first 
generations of stars will depend heavily on the reliabil- 
ity of these models. As these models mature and are 
(hopefully) shown to be reliable predictors of s-process 
nucleosynthesis at low metallicity, they will provide in- 
sight into the structure and evolution of stars on the 
AGB that cannot be obtained from the classical ad hoc 
approach to understanding the s-process. 

9.4. The Role of Future Studies of Neutron- Capture 

Isotopic Fractions 

We suggest that the features which are strongly 
blended with the n-capture lines in our stars (in par- 
ticular, Ti, V, Fe, and Ni) may be somewhat diminished 
in very metal-poor, n-capture enhanced stars. The re- 
quirements of high spectral resolution and high S/N to 
clearly discern n-capture isotopic fractions are currently 
difficult to achieve for the numbers of such stars that 
have been detected by the HK Survey dBeers et al.lll985l 
1992 1 ) and the Hamb urg/ESQ Survey (jWisotzki et all 
2000l IChristliebl 120031 ). owing to the faintness of these 
stars. From a technical outlook, observations of this 
nature will become more feasible with high-resolution 
echelle spectrographs on the next generation of 20— 30 m 
telescopes. Knowledge of the n-capture isotopic compo- 
sition of these stars would complement the information 
gleaned from elemental abundances alone and further en- 
hance our understanding of stellar nucleosynthesis and 
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chemical evolution throughout the early history of our 
Galaxy. 

10. CONCLUSIONS 

We have successfully measured the isotopic fractions of 
Eu and Sm in HD 175305, and we have assessed the Nd 
isotopic fraction as well. We find that the Sm isotopic 
fraction is suggestive of an r-process origin. The Eu and 
Nd isotopic fractions are unable to distinguish between 
an r- or s-process origin. Along with measurements of el- 
emental abundances of n-capture species, our Sm isotopic 
fraction reinforces the assertion that HD 175305 has been 
enriched by r-process nucleosynthesis. Both the elemen- 
tal and isotopic abundance distributions in HD 175305 
join the growing preponderance of evidence that supports 
the hypothesis of a universal r-process mechanism for el- 
ements with Z > 56. 

We have measured the isotopic fraction of Sm in 
HD 196944, which suggests an s-process origin. The Nd 
isotopic fraction is unable to distinguish between the two 
processes. The Sm isotopic fraction complements previ- 
ous elemental abundance measurements and observed ra- 
dial velocity variations of this star to confirm that the en- 
richment pattern of HD 196944 is consistent with transfer 
of s-process material from a companion star in the AGB 
phase of stellar evolution. 

We suggest that measurements of the Sm isotopic frac- 
tion will be more advantageous than measurements of 
the Eu isotopic fraction alone. The Sm isotopic frac- 
tion /152+154 permits a clearer distinction of pure-s- and 
pure-r-process content than Eu isotopic fractions can. 
There often are fewer blending features in the spectral 
regions where the useful Sm lines arc located, and there 
are more Sm lines available for analysis. Nd lines may 
offer more opportunities to measure an isotopic frac- 
tion in s-process-enriched stars than either Eu or Sm, 
although it is extremely challenging to discern the Nd 
isotopic fraction. We propose that measurements of 
isotopic fractions of multiple n-capture species in the 
same metal-poor star can allow direct quantification of 
the relative s- and r-process contributions to n-capture 
material, enable comparison of these contributions as a 



function of increasing nuclear mass number, test yields 
of low-metallicity stars on the AGB, and provide tests 
of s-process nucleosynthesis predictions. We have con- 
ducted one such test for 151 Eu using Eu and Pb abun- 
dances collected from the literature, finding that the 
lArlandini et alJ (|1999f ) stellar model predictions are a 
better match than the classical method predictions for 
151 Eu, although this conclusion hinges on measurements 
of the Eu isotopic fraction in only two (r + s)-enriched 
stars. 

Based on our results, we argue that the rare earth s- 
and r-process abundance patterns, which are observed 
in metal-poor stars at the elemental abundance level, 
are also present in the isotopic fractions. This implies 
that our understanding of n-capture nucleosynthesis is 
not wildly mistaken. This result has often been tac- 
itly assumed in studies of nucleosynthesis and chemical 
abundance patterns in the Galactic halo, and we contend 
that increasing empirical evidence for multiple rare earth 
species now exists to support this hypothesis. 
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APPENDIX 



COMMENTS ON INDIVIDUAL TRANSITIONS 

Here we describe the individual characteristics of our lines of interest as well as blending features. Where possible, we 
have tried to incorporate the highest quality atomic data for these blending features, including log(<?/) values calculated 
from laser-induced fluorescence measurements of the radiative level lifetimes and branching fractions measured from 
Fourier transform spectrometry. We alert the reader that many of the studies published before the early 1980s used 
less precise methods, and we employ those values with some caution. 



Europium Lines 

The Eu 11 line at 4129A is a clean, relatively unblended line in HD 175305; this spectral region was not observed 
in our other two stars. Th e only significant b lend is a Fe 1 line at 4129.4611 17 A in the far blue wing; th e wavelength 
of this line was tak en fromlNave et alJ (|1994h and the log(<?/) value was taken from lFuhr fe Wiesel (|2005l ). who scaled 
the value found by iMav et al.l ( 1974f ). We note that these values are not critical to determining the isotopic fraction 



of Eu. 

The Eu 11 line at 4205A is strongly blended with CH (and, to a lesser extent, CN) features, two V II lines, and a Y 11 
line. We adopt a CH linclist from the Kurucz database and a CN linelist from the Plez database 18 . Precise laboratory 

17 For wavelengths of blended features, we quote the number of digits to the right of the decimal point to the precision given in the 
reference and the precision used in our syntheses. 

18 Available online: ftp://saphir.dstu.univ-montp2.fr/GRAAL/plez/CNdata 
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wavelengths and log(g/) measurements are not available for the V II line at 4205. 04A; we adopt the V n 4205. 084A 
\og(gf) value from Bicniont et al.l (119891). A log(qf) va lue for the Y II line at 4204. 694A can be calculated from the 
oscillator strength presented in lHannaford et alj 1l982f). Performing a synthesis of this region of the Solar spectrum 
using an interpolate d empirical iHolweger fe Mullerl (|1974l ) model atmosphere, we find rather poor agreement with the 
iKurucz et al.l (|1984l ) Solar flux spectrum, with the most significant discrepancies arising from the CH bands and the 
wavelength of the Y n line. One would hope that in a very metal-poor star with sub-solar C abundance (HD 122563) 
these effects would be greatly diminished; the situation is much less hopeful for a C-enhanced star (HD 196944). Yet 
given the extent and severity of the blends of this Eu line, as well as its weakness in HD 122563, we are reluctant to 
trust any of the isotope fractions derived from this line. This is unfortunate because only this Eu line is strong enough 
to permit measurement of the elemental abundance in our observed spectra in HD 122563 and HD 196944. 

The Eu II line at 4435 A is strongly blended with a Ca I line at 4435. 67A and mildly blended with a Ni I line at 
4435. 33A. In HD 122563 and HD 196944, this Eu line is only a small sliver on the blue wing of the Ca line, and 
therefore no reasonable assessment of the elemental abundance or isotopic fraction can be made. In HD 175305 the Ca 
line is saturated, and we cannot produce a satisfactory fit to the blended profile no matter what log(<?/) value we adopt 
for this line. We adjust the log(gf) value of the Ni line to match the Solar spectrum, yet an additional adjustment of 
+0.65 dex to the Solar log(gf) value is necessary to bring this line into rational agreement with our observed spectrum. 
Even then we are not satisfied with the fit of the blue wing of the Eu line. Consequently, because such a small region of 
the spectrum covered by the Eu line remains unblended, we also disregard this line in our measurement of the isotopic 
fraction. 



Samarium Lines 

The Sm 11 line at 4424.34A is relatively strong and is ble nded only slightly with a Cr 1 line at 4424. 28A. An 
unpublished log(g/) value for this line, —0.63, was measured bv lSobeck et al.l (f2007f l . We measure Sm isotopic fractions 
from this line in HD 175305 and HD 196944; the line is too weak to measure well in HD 122563. 

The Sm 11 line at 4467 . 34A is also relativ ely strong and blended slightly with a weak Fe 1 line at 4467. 43A and an 
unidentified (in, e.g., the lMoore et al.lll966l atlas) line at 4467. 21 A. No experimental log(gf) value is available for the 
Fe line, but we use an inverted Solar analysis to derive a \og(gf) of —2.92 for this line. The unidentified feature can 
be reasonably modeled as an Fe line. The isotope shifts of this line are almost negligible and the HFS is also small, 
but lEkeland fe Hau gc (1975) used it in their analysis of the Sm isotopic fraction in the Sun This line is not covered in 
any of our stellar spectra, but we do use it to revisit the Sm isotopic fraction in the Solar photosphere. 

The Sm II 4591. 8lA line lies near the edge of the blue wing of a strong Cr 1 line at 4592. 05A that is easily accounted 
for in the synthesis. Of greater concern, however, is an apparently unidentified blend at 4591. 73A in HD 175305. We 
do not measure an isotopic fraction from this line but suggest that it might be useful in future studies. In HD 122563 
this line just falls off an echelle order, and in HD 196944 this line is too weak to observe. 

The Sm 11 4593. 53A line is affected by a very weak Fe I line at 4593. 53A, whose \og(gf) value is given bv lFuhr fc Wiesd 
(2005), which has been scaled from iMav et alJ (|1974h . We also notice a small, unidentified blend in the blue wing of 
this line at 4593. 42A; a close examination of iLundqvist et aTl (|2007h 's spectrum of CS 31082-001 in Figure 5 shows 
a similar feature. This Sm line is too weak to be observed in HD 122563, is observed but is too weak to assess 
the isotopic fraction in HD 196944, and, regrettably, in HD 175305 this line fell on a piece of the CCD where we 
encountered unrecoverable flat-fielding errors during the reduction process. 

The Sm 11 4595. 28A line is severe ly blended with a stro ng Fe I line at 4595.359lA. The \og(gf) value for this Fe line 
can be computed from the data in lO'Brian et al.l (|1991l ). but we are hesitant to believe any Sm abundances derived 
from this line in any of our stars. 

The Sm II line at 4604. 17A is promising. A Fe 1 line at 4603. 949A affects only the bluemost wing of this line, 
and a weak Fe 1 line blends with it at 4604. 24A; a log(gf) value is available for the latter of these two Fe lines in 
iGurtovenko fc Kostikl (|1981f ) and the log(gf) for the form er is derived from an inverted Solar analysis (—3.15), which 
provides a far better fit than the log(g/) rated an "E" in lFuhr fc W iese (20051) . We measure an isotopic fraction for 
this line in HD 175305, but the line is too weak in the other two stars. 

The Sm II 469 3. 63 A line is blende d with a Ti 1 line at 4693. 67A. Using our Ti 1 abundance and the log(gf) value 
for this line from lKuhne et a l. (1978), we cannot produce a satisfactory fit for the Sm line. An additional unidentified 
blend occurs at 4693. 75A in our spectrum of HD 175305. Our spectra of HD 122563 and HD 196944 do not cover 
this region. We do not make a measurement of the elemental abundance or isotopic fraction of Sm from this line in 
H P 175305, but w e sugg est that this line could be a candidate for future Sm isotopic analyses. 

ILundqvist et alJ (120071) note that the Sm 11 line at 4715. 27A is blended with a Ti 1 line at 4715.230A. A log(g/) 
value is available in lKuhne et al.l (|1978[ ). This line was not covered in any of our spectra, but isotopic analysis of this 
Sm line may be possible. 

The Sm 11 line at 4719. 84A is covered in all three of our spectra, but is only strong enough to measure in HD 175305 
and HD 196944. Full sets of hypcrfin e A and B constants ar e not available for this line, so we compute the profile using 
only the isotope shifts measured by Lundqvi st et alJ (120071) ■ This l ine is blended with a weak La 11 line at 4719. 92A; 
this La line was not included in the analysis of lLawler et al.l (|2001al ). While the Sm line is much stronger than the La 
line in both of our spectra, the uncertainty in the La log(g/) value and the lack of HFS structure information for the 
Sm line lead us to disregard this line in our isotopic analysis. 

The Sm 11 line at 5052. 75 A is not covered in any of our spectra. Full sets of hyperfine A and B constants are not 
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available for this line, either, but we compute a synthesis using the isotope shifts only. This line is blende d with a 
Ti I li ne at 5052. 87A and a Fe I line at 5052.9814A. A log(q/) value i s listed for the Ti line in the catalog of ISavanovl 
(1990), and a wavelength for the Fe line is given in lNave et alj (|1994h . We suggest that these blends may be able to 
be accounted for even without precise atomic data, and the isotope shifts alone could provide some indication of the 
Sm isotopic fraction in metal-poor stars, but measurement of an exact isotopic fraction is probably unwarranted here. 

The weak Sm n line at 5069. 47A is blended with a weak Fe I line at 5069.423 3A, which is blended with a weak Ti I 
line at 5069. 35A. The wavelength for this Fe line is listed in iNave et al.l (|1994f ). but we cannot find an experimental 
log(g/) value for this line, and the Ti blend prevents us from using an inverted Solar analysis to derive a log(gf) value 
for either line. We observe the Sm line in our spectra of HD 122563 and HD 196944; however, in both cases the Sm 
line is weak and the observed spectrum is littered with small, undocumented blends in addition to the Fe line noted 
here. We recommend that any attempts to use this line for Sm isotopic analysis proceed with great caution. 

The Sm II line at 51 03. 09A is weak and li es on the red edge of a strong Ni I line at 5102. 964A. A log(g/) value for 
the Ni line is given in iDoerr fc Kocl (fl985h . and we have set the wavelength of this line from analysis of the Solar 
spectrum. Full sets of hyperfine A and B constants are not available for this line, yet we compute a synthetic spectrum 
using only the isotope shifts. We can detect the Sm line in HD 175305, but it is so overwhelmed by the Ni line that we 
cannot trust even our derived elemental abundance. This Sm line is too weak to det ect in HD 122563 and HD 196944. 

Our redmost Sm II line lies at 5104. 48A. It is blended by a Fe I line at 5104.43 75A. lNave et all (fl994l) quote a precise 
wavelength for this Fe line, and its log(gf) value is listed in lFuhr fe Wiesd (|2005D , who scaled lMav et al.l (|1974l )'s value. 
When we synthesize this line in HD 175305, we are forced to enhance the Sm abundance by more than a factor of 2. 
We have re-evaluated th e logf q f ) value of th e Fe blend by comparing with the Solar spectrum, and find that we can 
reasonably only increase iFuhr fe Wiesd (|2005f )'s \og(gf) value by « 0.2 dex, which does little to provide a better match 
between our observed and synthetic spectra. In addition, there is a significant unidentified blend at 5104. 59A in the 
red wing of our Sm line. While it is clear that the Sm line is broadened by HFS and isotopic shifts, we unfortunately 
are unable to measure an isotopic fraction. This Sm line is too weak to detect in HD 122563 and is not covered in our 
spectrum of HD 196944. 



Neodymium Lines 

The Nd n line at 4177. 32A is relatively strong, and only its far red wing is blended w ith a strong Y II and Fe I blend. 
A log(gf) value can be calculated for the Y II line at 4177. 536A from the /-value in lHannaford et ail (|1982h . and a 
log(g/) value for the Fe I line at 4177.5935A is given in lO'Brian et al.l (|1991h . This blend is easily accounted for and 
has little effect on the derived Nd isotopic fraction. This Nd line was not covered in our spectra of HD 122563 and 
HD 196944, but we measure an isotopic fraction from this line in HD 175305. 

The Nd n line at 4232. 38A has two identified blends. The Hf II line at 4232. 386A lies coincident with our Nd 
line; a log(g/) value for this line was measured by iLundqvist et"aT (|2006h . A V I line at 4232. 46A also blends with 



our Nd line, and a log(gf) value is provided by iMartin et alJ ( 1988 ) in the NIST database. The uncertainties in the 



abundances and transition probabilities of these two lines, as well as their severity of blending with our Nd feature, 
lead us to disregard this line in our spectrum of HD 175305. This line was not covered in our spectra of HD 122563 
and HD 196944. The isotope shifts of this Nd line are noticeable, though, so we report the hyperfine patterns for this 
line nonetheless. 

The Nd n line at 4314.5lA lies in a region where continuum placement is difficult due to the presence of the CH 
G-band. We have not attempted to measure the Nd isotopic fraction from this line in any of our stars, although the 
Nd line is visible in HD 175305 and HD 196944. The isotope shifts of this line are rather large (~ 0.035A), so we 
report the hyp erfine patterns for t his lin e, in hopes that the Nd isotopic fraction could be discerned from this line in a 
favorable star. iDen Hartog et all (|2003f ) did not measure a log(gf) value for this line, but we derive \og(gf) — —0.22 
from an inverted Solar analysis, assuming loge (Nd) Q = 1.50. 

The Nd II line at 4358. 16A is relatively strong and unblended in HD 175305, and we use it to measure an isotopic 
fraction. This line was not covered in our spectra of HD 122563 and HD 196944. 

The Nd II line at 4446. 38A has the largest isotope shifts of any of the Nd transitions for which we report an 
isotopic fraction. I n HD 175305, this line is blended only with a weak Gd II line at 4446. 502A in the red wing, and 
IDen Hartog et all (|2006l ) provide a log(gf) value for this line. In HD 196944, some unidentified blends are observed 
in both the red and blue wings of this line, but it does not appear that these blends affect our measurement of the 
isotope fraction. We measure an isotopic fraction from this line in each of these two stars. This line is observed in our 
high-S/N spectrum of HD 122563, but it has a continuum depth of only ~ 1.5%, and we cannot measure an isotopic 
fraction from this line. 

The Nd II line at 4567. 6lA has a depth of only 4% in HD 175305, and an isotopic fraction can only marginally be 
deduced from this line. We identify no blending features, although our observed spectrum clearly shows some weaker 
blends in the wings of this line. This line was not covered in our observed spectra of HD 122563 and HD 196944. 

HYPERFINE COMPONENT DATA FOR SAMARIUM 

We present hyperfine and isotopic components for the 13 lines of Sm n described in § IA.2l in Table |4] Hyperfine A 
(magnetic dipole) and B (electric quadrupole) constan ts, which govern the rela tive positions of the hyperfine compo- 
nents, along with isot ope shifts, were taken first from iMasterman et all (|2003[) if avail a ble. C onstants for additional 
levels were taken from iDorschel et ail ()198lD . lYoung et alj (|1987l ), and IVillemoes et alj ([1995). We give preference to 
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the hyper fine constants and i sotop e shi fts determined fr o m rad io frequency measurements or with a single- frequency 
laser, as i nlMasterman et al.l (120031 ) and iBeiersdorf et al.l (|1995T ) . supplemented with isotope shifts measured from FTS 
spectra in lLundqvist et al.1 (|2007l ) for lines not covered in earlier studies. 

The relative strengths are calculated from the LS (Russell-Saunders) angular momentum coupling formulae (e.g., 
iCondon fc Shortlevlll95l p. 238), where we have replaced the orbital angular momentum quantum number (L) by 
the total electronic angular momentum (J), the electron spin (S 1 ) by the nuclear spin (I), and the total electronic 
angular momentum (J) by the total atomic angular momentum (F). Thus, F = I + J, and the energy of the J 
level is split into a number of sublevels given by F, which runs from I — J\ to / + J. The energy spacing between 
each component is proportional to F. The nucleus may also have an electric quadrupolc moment, which produces 
an additional shift (but not splitting) in the energy levels. The strongest relative components are those for which 
A F = A J. We normalize the strengths for each line such that they sum to one. Absolute transition wavenumbers 
are from iLunda vist et~l . (2007) i f the l ines were included in their study, otherwise they are computed using the energy 
levels tabulated bv lMartin et al.l (|1978l ). The center-of-gravity transition wavenumbers and air wavelengths are given 
in the first two column s , resp ectively. Component wavenumbers are converted to air wavelengths according to the 
formula given bv lEdlir] (fl953h . 

HYPERFINE COMPONENT DATA FOR NEODYMIUM 

We present hyperfine and isotopic components f or the 6 lines o f Nd II described in § IA.3I in Table [5] Hyperfine A 
and B constants and isotope shifts were taken from lRosner et all (|2005f ). All calculations were performed analogously 
to the Sm case. Cente r-of-gravity transition wavenumbers for the lines of interest were taken from FTS measurements 
bv lBlaise et "all (|1984f ). 
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TABLE 1 

Comparison of Model Atmosphere Parameters 



Reference 


^ CII 


loff CI ) 


fM/Hl 


[Fc/Hl a 


• 

^micro 


"macro 












(km s ) 


(km s ) 






HD 


122563 








This study 


4430 


0.55 


-2.63 


-2.83 


2.4 


4.4 


Fulbright (2000) 


4425 


0.6 




-2.6 


2.75 




Johnson (2002) 


4450 


0.50 


-2.65 


-2.75 


2.30 




Honda et al"T2004b) 


4570 


1.1 




-2.77 


2.2 




Simmcrcr et al. ( 2004) 


4572 


1.36 




-2.72 


2.90 








HD 


175305 








This study 


4870 


2.15 


-1.40 


-1.60 


1.2 


4.6 


Fulbrieht (2000") 


4575 


2.4 




-1.30 


1.4 




Burris et al. (2000) 


5100 


2.50 


-1.50 


-1.40 


1.40 




Cowan et al. (2005) 


5040 


2.85 




-1.48 


2.00 








HD 


196944 








This study 


5170 


1.80 


-2.26 


-2.46 


1.7 


6.8 


Zacs et al. (1998) 


5250 


1.7 


-2.0 


-2.45 


1.9 




Aoki et al. (2TO2) 


5250 


1.8 




-2.25 


1.7 





a [Fe i/H] is listed here. We find [Fe Il/H] = -2.79 for HD 122563, 
[Fe n/H] = -1.62 for HD 175305, and [Fe o/H] = -2.43 for HD 196944. 
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TABLE 2 

Isotope Mixes for Pure s- and t--process Material and Measured Isotopic Fractions 

stellar model a classical method 
Species s-only r-only s-only r-only Solar System HD 175305 HD 196944 



Ba: / 0( jd = 


0.11 


0.46 


0.09 


0.72 


0.178 






Nd: /142+144 = 


0.67 


0.32 


0.69 


0.27 


0.510 


21+ - 56 


S6+ 64 


Sm: /152+154 = 


0.21 


0.64 


0.20 


0.64 


0.495 


0.51 ±0.08 


0.35 ±0.14 


Eu: /i 5 i = 


0.54 d 


0.47 


0.04 d 


0.47 


0.478 


0.50 ±0.04 





a IArlandini et al.l (fl999h b ISlmmerer et all l[200l ) : [Cowan et al.l p00^ c IRosman fc Tavlorl (fl998l ) d The 
classical fit — based on assuming a smoothly-varying abundance cu rve^ — predicts a much smaller 
s-process contribution to 151 Eu than does the stellar s-process model of Arlandini et al. (1999). See §[8] 
for further discussion. 
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TABLE 3 

Heavy ti-capture Elemental Abundances in HD 175305 



l°g(9/) lo S e© 



Species 


Z 


log e 


[X/Fe] 


a 


1 ' lines 


log £@ 


Ref. 


Ref. 


Ba II 


56 


+0.91 


+0.40 


0.16 


1 


2.13 


1 


10 


La II 


57 


-0.29 


+0.20 


0.13 


3 a 


1.13 


2 


2 


Ce II 


58 


+0.05 


+0.09 


0.13 


3 


1.58 


3 


10 


Nd II 


60 


+0.09 


+0.21 


0.13 


10 


1.50 


1 


10 


Sm II 


62 


-0.25 


+0.37 


0.11 


9 


1.00 


5 


5 


Eu II 


63 


-0.65 


+0.46 


0.07 


3 


0.51 


6 


10 


Gd II 


64 


+0.00 


+0.50 


0.15 


4 


1.12 


7 


10 


Dy II 


6(5 


+0.01 


+0.49 


0.16 


1 


1.14 


8 


10 


Pb I 


82 


< +0.10 


< +0.45 




1 


1.95 


9 


10 



References. (1) IKlose et al.l [I2002T); ( 2) ILawler et all (I2001a|) ; 

(3l lPalmeri et al l II2000I) ; (4 ) IDen Hartoe et al] d2003l> ; ( 5) ILawler et al] (12006 ); 
(6) ILawler et al.l 12001blT ; (71 IDen Hartoe et al] (120061); (8 ) [Wickliffe et~aLl (125000 ; 
O nBiemont et al.l (120001 ); f 10) IGrevesse fc Sauvaf(l2002T) 

Note. — Ionized species abundance ratios are referenced to the Fe II abun- 
dance in Table [3] We adopt t he Solar photospheric Fe abundance logeQ (Fe) = 
7.52 from lSneden et al] l tl99lft . 

a We synt hesize these La II lines using the HFS patterns listed in Table 4 of 
llvans et al] ((2006] '). 
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TABLE 4 

Isotopic Hyperfine Structure Patterns for Samarium 



Wavenumbcr 
(cm -1 ) 



Wavelength Component Position 
(A) (cm" 1 ) 



Component Position Isotope 
(A) 



Strength a 



22595.910 
22595.910 
22595.910 
22595.910 
22595.910 
22595.910 



4424.3373 
4424.3373 
4424.3373 
4424.3373 
4424.3373 
4424.3373 



0.20775 
0.13532 
0.06911 
0.15704 
0.09083 
0.03117 



-0.04067 
-0.02649 
-0.01353 
-0.03074 
-0.01778 
-0.00610 



117 
147 
147 
117 
147 
147 



0.028106 
0.001521 
0.000040 
0.023908 
0.002540 
0.000097 



Note. — The full table is available in machine-readable form in the electronic edition of the 
journal; only a small portion is shown here to present its general form and content. For lines 
without full sets of hyperfine A and B constants available (4719, 5052, and 5103), for the 147 Sm 
and 149 Sm isotopes we list only the isotope shifts, computed as if all A and B values are zero. The 
isotope shift for t he 14 4 Sm isotope of the 5069 line could not be measured in the FTS spectrum of 
Lundavist et al. (2007). We estimate this value by scaling the energy shifts of the 144 Sm isotope 
relative to the heavier isotopes from the 5052 and 5103 lines, which show similar structure patterns 
(M. Lundqvist, 2007, private communication). 

a The sum of the strength values for each line have been normalized to one. When using these 
components in ou r spect ral synthesis, we divide out the Solar System isotopic ratios as given by 
Rosman & Tavlor (1998), but we stress that this step is NOT included in the data presented here. 
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TABLE 5 

Isotopic Hyperfine Structure Patterns for Neodymium 



Wavenumbcr 
(cm" 1 ) 



Wavelength Component Position Component Position Isotope 
(A) (cm" 1 ) (A) 



Strength a 



23932.050 
23932.050 
23932.050 
23932.050 
23932.050 
23932.050 



4177.3196 
4177.3196 
4177.3196 
4177.3196 
4177.3196 
4177.3196 



0.10001 
0.14184 
0.05369 
0.17828 
0.09013 
0.01208 



-0.01746 
-0.02476 
-0.00937 
-0.03112 
-0.01573 
-0.00211 



143 
143 
143 
143 
113 
143 



0.024146 
0.001719 
0.019886 
0.000061 
0.002841 
0.016161 



Note. — The full table is available in machine-readable form in the electronic edition of the 
journal; only a small portion is shown here to present its general form and content. 
a The sum of the strength values for each line have been normalized to one. When using these 
components in ou r spect ral synthesis, we divide out the Solar System isotopic ratios as given by 
Rosman fc Tavlorl ^19981 1 . but we stress that this step is NOT included in the data presented here. 
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1 44 

3m 

143 
Pm 



0% 



145 

Sm 

144 
Pin 



14C) 
Sm 

145 
Pm 



147 

Sm 

146 
Pm 



92% 



148 

Sm 

147 

Pm 



0% 



149 
Sm 

148 
Pm 



86% 



151 
lu 
I 

150 

Sm 

149 
Pm 



100% 



0% 
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Fig. 1. — A small region of the table of nuclides, showing increasing mass number (A) on the horizontal axis and increasing atomic 
number (Z) on the vertical axis. Naturally-occurring isotopes are represented by the gray squares. The figure highlights the naturally- 
occurring isotopes of neodymium (Nd, Z = 60), samarium (Sm, Z = 62), and europium (Eu, Z = 63); promethium (Pm, Z = 61) has no 
na turally-occurr i ng iso topes. We also show the predicted percentage of r-process material in the S. S. abundance for each isotope as given 
by Cowan ct al. (2006). This figure illustrates the pathways by which the s- and r-processes operate and the resulting differences in the s 
or r fraction of S. S. material. 
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Fig. 2. — One method used to match the observed and synthetic spectra, displaying how the point of insensitivity has been matched 
to the observed spectrum. The pure-s- and pure-r-syntheses are indicated by dot-dashed and dashed lines, respectively. The observed 
spectrum is represented by the histogram. The point at which the pure-s- and pure-r-process syntheses intersect is defined to be the point 
of insensitivity to the isotopic mix. We adjust the wavelength of the synthetic spectra so that the point of insensitivity coincides with the 
center of the nearest pixel in the observed spectrum. The isotopic fraction can then be measured using the fitting algorithm. The vertical 
sticks represent the positions and relative strengths of the individual hyperfine components; solid sticks represent the 152 Sm and 154 Sm 
isotopes, while dashed sticks represent the remaining five isotopes. 
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Fig. 3. — Heavy n-capture elemental abundan ces for HD 17530 5. Black squares indicate our observed abundances. Gray triangles 
indicate the o bserv ed Ba and Eu abundances from [Fulbright (2000). Gray circles indicate the observe d Os, Ir, and Pt abundances from 
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Fig. 4. — Our synthesis of the Eu II line at 4129A in HD 175305. We have measured an isotopic fraction of /151 = 0.50 ± 0.04 for this 
line. The observed spectrum is indicated by open squares. In panel a, we show the syntheses with the extreme values of /151 as well as the 
pure-s- and pure-r-process syntheses. In panel b, the solid curve represents our best fit isotopic fraction, while the dotted curves represent 
the 3cr uncertainties. In both panels, the long-dashed curve is a synthesis with no Eu present, indicating that this line is generally free of 
blending features. The vertical sticks represent the positions and relative strengths of the individual hyperfinc components; dotted sticks 
represent 151 Eu while solid sticks represent 153 Eu. Note that the Eu elemental abundance reported is the best-fit value for this synthesis 
and not our derived Eu abundance for HD 175305. 
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Fig. 5. — Ax 2 values for the Eu 4129 line in HD 175305, shown for three values of Vmacro- The solid line is a third-order polynomial 
fit to the measured A% 2 values for the best value of Vmacro, 4.6kms _1 , and the dotted lines are fits to the measured Ax 2 values for the 
uncertainties in Vmacro- One, two, and three-cr confidence intervals are indicated by the dashed lines. The fits have well-defined minima, 
and the choice of the macroturbulent velocity has little effect on the derived isotopic fraction. 
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Fig. 6. — Our synthesis of the Eu II line at 4205A in HD 175305. Symbols are the same as in Figure[4] We emphasize that the isotopic 
fraction shown in panel b is the value derived from the 4129A line and not a fit to this line. Note the strong blends. The isotope fraction 
derived from the 4129A line provides a reasonable fit to the 4205A line as well, requiring only small changes in the elemental abundance. 
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Fig. 7. — Our synthesis of the Eu II line at 4435A in HD 175305. Symbols are the same as in Figure[4] We emphasize that the isotopic 
fraction shown in panel b is the value derived from the 4129A line and not a fit to this line. Note the strong blend with the saturated Ca I 
line at 4535. 67A. The isotope fraction derived from the 4129A line provides a reasonable fit to the 4435A line as well, requiring only small 
changes in the elemental abundance. 
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Fig. 8. — Our synthesis of the Sm II line at 4424A in HD 175305. The observed spectrum is indicated by open squares. In panel a, 
we show an (O — C) plot for the pure-s- and pure-r-process syntheses shown in panel b. The unshaded region indicates the wavelengths 
over which we measure the isotopic fraction using the x 2 algorithm. In panel b, we show syntheses for the pure-s- and pure-r-process 
syntheses (the dot-dashed and short-dashed curves, respectively) illustrating the contrast between these extremes. In panel c, we show the 
synthesis of our best fit, represented by the solid curve, and the 3<r uncertainties, represented by the dotted curves. The long-dashed curve 
is a synthesis with no Sm present, indicating the presence of blending features. The vertical sticks represent the positions and relative 
strengths of the individual hyperfinc components; solid sticks represent the two heaviest isotopes, while dashed sticks represent the lighter 
five isotopes. We measure a Sm isotopic fraction of /152+154 = 0.551q'^3 for this line, which suggests an r-process origin. Note that the 
Sm elemental abundance reported is the best-fit value for this line only and not our mean [Sm/Fe] in HD 175305. 
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Fig. 9. — Our synthesis of the Sm II line at 4424A in HD 196944. Symbols are the same as in Figure [8] We measure a Sm isotopic 



fraction of /152+154 



0.35 



+0.11 



for this line, which moderately suggests an s-process origin. Note that the Sm elemental abundance 



reported is the best-fit value for this line only and not our mean [Sm/Fe] in HD 175305. 
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Fig. 10. — Our synthesis of the Sm II line at 4604A in HD 175305. Symbols are the same as in Figure [8] We measure a Sm isotopic 



fraction of /152+154 



0.50 



+0.09 



from this line, which suggests an r-process origin. Note that the Sm elemental abundance reported is the 



best-fit value for this line only and not our mean [Sm/Fe] in HD 175305. 
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Fig. 11. — Ax 2 values for the Sm 4604 line in HD 175305, shown for three values of V ma cro- Symbols are the same as in Figure IS1 
Here the macroturbulence has only a small effect on the derived isotopic fraction. Similar results are found for the Sm 4424 line in both 
HD 175305 and HD 196944. 
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Fig. 12. — A summary of our best isotopic fraction measurements for each species. The Nd isotopic fractions reported here were measured 
by the method of matching the observed and synthetic spectra at the point of insensitivity to the isotopic mix. The Sm and Eu isotopic 
fractions were measured by fitting the shape of the line profile. Squares represent measurements in HD 175305, while circles represent 
measurements in HD 196944. The dotted line and dashed line ( colored blue and red, respectively, in the electronic edition) represent the 
pure-s- and pure-r-process predictions of Arlandini ct al. (1999). The differences between the pure-s- and pure-r-process predictions are 
scaled together to emphasize the relative precision with which the two processes can be distinguished by each species. [See electronic edition 
for a color version of this figure.] 
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Fig. 13. — A summary of our best isotopic fraction measurements for each species. The data and symbols are the same as in Figure IT2l 
The measurements are cast here in an absolute sense, emphasizing the relative precision with which each species was measured. [See 
electronic edition for a color version of this figure.] 



8<S 



Roederer et al. 



1.0 - 



o.o - 



X) 

\ 

ia-1.0 

CO 

60 
O 



-2.0 - 



: ,i : 



-3.0 - 



T 



T 



O Aoki et al. (2001) 
A Aoki et al. (2002) 
X Aoki et al. (2006) 
□ Barbuy et al. (2005) 
O Cohen et al. (2003) 
O Cohen et al. (2006) 
# Cowan et al. (2002) 



■1 



# Frebel et al. (2007) -< Jonsell et al. (2006) 

▲ Hill et al. (2002), % Masseron et al. (2006) 

Aoki et al. (2003a). O Preston et al. (2006) " 

Plez et al. (2004) + Sneden et al. (2003) 

-X- Ivans et al. (2005) ■ Yong et al. (2006) 

+ Ivans et al. (2006) ▼ this study 
V Johnson & Bolte (2002) 

Arlandini et al. (1999) 

= Cowan et al. (2006) 



— pure-r 





-3.5 



-3.0 



-2.5 
[Fe/H] 



-2.0 



-1.5 



Fig. 14. — The log e(Eu/Pb) abundance as a function of [Fe/H] for 36 halo and globular cluster stars collected from the literature or 
measured in our study. References for the individual stars are indicated in the figure. Solid shapes indicate stars classified as r-process 
enriched, open shapes indicate stars classified as s-process enriched, and starred shapes indica te stars classifi ed as (r + s)-process enric hed. 
Upper limits are indicated with downward arrows. The abundance predictions from Arlandini et al. (1999) and Cowan et al. (2006) are 
shown as dotted and dashed lines, respectively. (In the electronic edition, blue and red lines represent the pure-s-process and pure-r-process 
abundance predictions.) The s- and (r + s)-process-enriched stars cannot unambiguously distinguish between the different pure-s abundance 
predictions. [See electronic edition for a color version of this figure.] 
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Fig. 15. — The log e( 151 Eu/Pb) abundance as a function of [Fe/H] for 7 halo stars collected from the literature or measured in our 
study. References for the individual stars are indicated in the figure. Symbols are the same as in Figure [14] The two (r + s)-enriched stars 
favor the pure-s 151 Eu abundance prediction of Arlandini et al. (199SJ. [See electronic edition for a color version of this figure.] 
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Fig. 16. — The log e (Ba/Nd) abundance as a function of [Fe/H] for 26 halo stars collected from the literature or measured in our study. 
References for the individual stars are indicated in the figure. Symbols are the same as in Figure 1141 No distinction is observed between 
the s- and (r + s)-enriched stars from their Ba and Nd abundances, leading us to conclude that the two (r + s)-enriched stars in Figure [TS] 
are representative of stars containing products of predominantly-s-process-enriched material. [See electronic edition for a color version of 
this figure.] 



